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Chapter 1 
Introduction 



Several possible scenarios of proton - nucleus reaction are considered 
nowadays. According to one of the scenarios, incoming proton deposits en- 
ergy into target nucleus. It knocks out a few nucleons and leaves excited 
residual nucleus. Then, nucleons and various fragments are emitted from 
the excited residuum. This scenario is called spallation. It is also possible, 
that the residuum splits up slowly (fissions) into two parts, which then emit 
particles. This scenario is known as emission from fission fragments. But, 
it could be also that all fragments appear simultaneously. This would have 
features of a phase transition in nuclear matter and is called fragmentation. 
These reaction scenarios are based on experimental observations of differ- 
ent final states. Generally, observation of one heavy nucleus (in respect to 
the mass of initial target), a small number of light fragments and numerous 
individual nucleons indicates spallation. Detection of a large number of in- 
termediate size fragments indicates fragmentation. Nevertheless, spallation 
and fragmentation are correlated. Their differentiation is not clear, difficult 
and still under discussions. Previous studies indicate, that spallation is the 
most probable scenario of proton - nucleus reaction, what will be shown also 
in this dissertation. The main idea of this work concerns theoretical study 
of proton induced spallation reactions in wide range of incident energy and 
mass of target nuclei; fission and fragmentation are not discussed in details. 

The following definition of spallation process can be found in Nuclear 
Physics Academic press: "Spallation - a type of nuclear reaction in which 
the high-energy of incident particles causes the nucleus to eject more than 
tree particles, thus changing both its mass number and its atomic number." 
So, the term spallation means a kind of nuclear reactions, where hadron with 
high kinetic energy (100 McV up to several GeV) interacts with a target. 
First, this term was connected with observation of residuum of reaction cor- 
responding to losses of mass of target nucleus from few up to several dozen 
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nucleons. Nowadays, it means mechanism, in which high energy hght parti- 
cle causes production of numerous secondary particles from target nucleus, 
leaving cold residuum of spallation. As a result of such process also various 
Intermediate Mass Fragments (IMF), i.e. fragments with masses in range 
4 < MjMF < 20, are observed. 

From historical point of view, the possibility of heating a nucleus via bom- 
barding by neutrons was suggested first time in 1936 by N. Bohr p. Studies 
of similar reactions were possible due to development of accelerator technics. 
It was in the end of fourties, when accelerators could provide projectiles with 
energies higher than 100 MeV [2]. Experimentally, two - component spec- 
tra of emitted particles are observed: anisotropic high energy part, which 
dominates in forward angles (i.e. the high energy tail decreases at backward 
angles, as it is seen in the example Fig. 11.11) and isotropic, low energy part. 
These general features of spallation process are established experimentally. 
A theoretical picture of an incident particle colliding successively with sev- 
eral nucleons inside target nucleus, losing a large fraction of its energy was 
proposed by Berber in 1947 [3]. Before, in 1937 Weisskopf considered possi- 
bility of emission of neutron from excited target nucleus [1]. In the end of 
fifties. Metropolis [5] and Dostrovsky [6] (who used the ideas of Serber and 
Weisskopf) suggested description of spallation as two step process involving 
energy deposition and subsequent evaporation. They formulated and per- 
formed first Monte Carlo calculations of the reactions. Such treatment of 
spallation reactions is used from that time up to now. In more details, the 
first, so-called fast stage of the spallation is highly non-equilibrated process. 
High energy proton causes an intra-nuclear cascade on a time scale ~ 10~^^s. 
The incident projectile goes through the target nucleus and deposits a sig- 
nificant amount of excitation energy and angular momentum, while ejecting 
only a few high energy nucleons and, with a minor yield, pions and light ions. 
The result of the first stage is excited residual nucleus in thermodynamical 
equilibrium (totally or partly equilibrated), with excitation energy ~ a few 
MeV/nucleon. 

In case of thick target, i.e. system of several nuclei, the ejectiles, as secondary 
projectiles can cause so-called inter-nuclear cascade, placing individual nuclei 
into excited states, as illustrated in Fig. 11.21 

The second, so-called slow stage of the spallation, consists in deexcitation 
of the residuum by evaporation of particles. The isotropic emission (in the 
system of nucleus) of nucleons (mainly neutrons), light and heavy ions (d, t. 
He, Li, Be, B, 7) takes place on a time scale ~ 10^^^ - lO^^^s. 

From many years spallation reactions of medium and high energy protons 
with atomic nuclei are still of interest for many reasons. First of all, because 
knowledge of the reaction mechanism is still not complete. This is interesting 
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Figure 1.1: Experimental observations of an example 2.5 GeV proton induced 
spallation reaction on Au target; inclusive alpha particles spectra measured 
by PISA at COSY 0? 



both from theoretical and experimental point of view. Experimental data of 
double differential cross sections of emitted particles in the reactions are nec- 
essary for testing, validation and developing of theoretical models. It means, 
experimentally measured cross sections for exclusive elementary reactions 
(e.g. NN, Nvr, ...) are implemented in theoretical models. Then, results 
of calculations are compared with results of inclusive measurements. It is 
reasonable to study the reaction mechanism on the base of proton - nucleus 
rather than nucleus - nucleus collisions, where all processes start to be much 
more complicated (e.g. presence of distortions due to collective processes 
like compression, deformation, high spin [8]). Moreover, proton - nucleus 
reactions are important and indispensable also for experiments of nucleus - 
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nucleus collisions (e.g. HADES [9], CHIMERA [lO]). Results of proton - 
nucleus reactions facilitate extraction and interpretation of results of nucleus 
- nucleus reactions. Other reasons concern very broad range of applications 
(e.g. in medicine (radiation therapy), cosmology, accelerator technology). 
Relatively huge number of produced neutrons suggested the idea of using 
spallation reactions as neutron sources. Nowadays, neutron beams are pro- 
duced in nuclear reactors. Reactors dedicated for such production generate 
also a lot of heat; about 190 MeV of energy is dissipated for single produced 
neutron. In accelerator based sources, neutrons are produced in a spalla- 
tion process, with only about 30 MeV of energy dissipated for one generated 
neutron. During last decade several spallation sources (IPNS [llj, ISIS [IT] . 
LANSCE [12], SINQ [13]) became operational. 

Spallation reactions are very important in accelerator technology (e.g. acti- 
vation of detectors, radiation protection). The reactions are used for energy 
amplification, also for production of energy from nuclear waste and further- 
more, transmutation of long - lived radioactive nuclei of nuclear waste to 
stable or short - lived, in order to avoid their long term storing |14j . 
Astrophysical models have to include spallation processes. If one compares 
abundances of cosmic rays and solar system elements, it is seen that Li, Be 
and B in cosmic rays are enriched by more than 6 orders of magnitude, as 
shown in Fig. 11.31 They were evidently produced in spallation reactions of 
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Figure 1.3: The relative abundances of elements in cosmic rays measured at 
Earth and in the solar system, taken from ji6| 



hydrogen nuclei (which consist about 87 % of cosmic rays) with heavy ele- 
ments (produced due to stars explosions). For more informations see p!5llT6]. 

Theoretical predictions of the process are important in each of mentioned 
above cases. Several models have been constructed in order to describe the 
spallation process. First stage of the reaction is described by a class of 
microscopic models, e.g. [171 [TSl HH] • For the second stage statistical models 
are used, e.g. [201 EI] ■ Nevertheless, the theoretical investigations are rather 
fragmentary and still not satisfactory. Global qualitative and quantitative 
informations of the spallation reactions are needed. 

Description of global average properties of proton induced spallation re- 
actions in wide range of projectile energy and mass of target nuclei is the 
main subject of this work. This is investigated within a transport model 
based on Boltzmann - Uehling - Uhlenbeck (BUU) equation supplemented 
by a statistical evaporation model. The used transport model has been spe- 
cially developed in the frame of this work, in order to enable description of 
such reactions. 

This work is organized as follows. In Chapter 2, present knowledge of the 
reaction dynamics and review of various theoretical models is given. Chap- 
ter 3 includes description of the Hadron String Dynamics (HSD) approach 
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of the first stage of the spallation reaction. In Chapter 4 and Chapter 5, 
results concerning properties of the first stage of the reaction are presented 
and discussed. Chapter 6 concerns pions, which could be produced only in 
the first stage of reaction. In Chapter 7, statistical evaporation models for 
the second stage are recalled. Chapter 8 contains bulk model predictions of 
both stages of proton - nucleus reactions. In Chapter 9 and Chapter 10, the 
models results and comparison with available experimental data and with 
results of other models are presented, respectively. Finally, in Chapter 11, a 
summary and conclusions are given. 



Chapter 2 



Present knowledge of the 
reaction dynamics - basic 
theoretical models 

Our understanding of physical phenomena is expressed as modelhng. At 
present, the broadest platform for such modelling is quantum mechanics ap- 
proach. Unfortunately, many - body systems are usually an extreme challenge 
for existing methods of quantum mechanics. One has to rely on rather sim- 
ple, much more straight formed concepts. In this Chapter, such general basic 
concepts will be shown as ingredients of typical models of nuclear reaction. 

Several microscopic models have been constructed in order to describe 
the first stage of proton - nucleus reaction. All of them have the same basis, 
they describe the reaction as a cascade of nucleon - nucleon collisions, but 
employing different assumptions. The main difference concerns implemented 
potential of nucleon - nucleus interaction. One can distinguish the simplest 
models, which neglect features of the mean field dynamics and employ con- 
stant static potential, like a class of Intra - Nuclear Cascade (INC) mod- 
els. Other, more sophisticated approaches comprise dynamically changing 
field and minimal fluctuations obtained due to use of test particle method, 
i.e. models based on Boltzmann - Uehling - Uhlenbeck (BUU) transport 
equation. There are also models, which include real fluctuations and parti- 
cles correlations, employing two- and three- body potentials, e.g. Quantum 
Molecular Dynamics models. The main ideas of the different types of nuclear 
reaction approaches are described below. 

In this Chapter, devoted for model-like approach of investigation of nu- 
cleon - nucleus reactions, also so-called percolation model is presented. This 
rather basic model describes fragment mass distributions very well (Section 

ESD- 



9 



10 



2.1 Intranuclear Cascade model 

Intranuclear cascade (INC) model is used as a base in many existing codes 
for the first stage calculations. Description of the nucleon - nucleus reaction 
in terms of binary nucleon - nucleon collisions inside nucleus is the basic as- 
sumption of the model. In principle, the single particles approach of the INC 
is justified as long as the de Broglie wavelength A of the cascade particles 
is smaller than the average internucleon distance in the nucleus (~ 1.3 fm). 
This indicates the low energy limit of the model (e.g. for projectile kinetic 
energy Tp = 100 MeV ^ A = 2.7 fm, for Tp = 1000 MeV ^ A = 0.7 fm). 
The INC calculations follow the history of individual nucleons that becomes 
involved in the nucleon - nucleon collisions in a semi - classical manner. 
It means, the momenta and coordinates (trajectories) of the particles are 
treated classically. The only quantum mechanical concept incorporated in 
the model is the Pauli principle. 

The first code of INC has been created by Bertini [22], in 1963. Later, the 
conception was used also in other codes, e.g. by Yariv in his ISABEL code 
[23] • In the 80's and 90's, the next versions of INC model was developed by 
Cugnon et al. [T7j . 

The main features of the standard INC approach are the following. 

The initial positions of target nucleons are chosen randomly in a sphere 

1 /3 

of radius R = 1.12 ■ Arp fm, where: At is the mass number of target nu- 
cleus. Momenta of the nucleons are generated inside a Fermi sphere of radius 
Pf = 270 MeV/c. Neutrons and protons are distinguished according to their 
isospin. All nucleons are positioned in a fixed and constant, attractive poten- 
tial well of Vo = 40 MeV depth, inside the nuclear target volume. The depth 
value is taken a bit higher than the Fermi energy {Ep ^ 38 MeV), so that 
the target is stable during the reaction. The idea of fixed average potential is 
based on a relatively low number of particles emitted during the INC stage of 
reaction, which disturbes the mean field only slightly. The incident particle 
(of incident energy Tp) is provided with an impact parameter b, chosen ran- 
domly on a disc of radius R. It is positioned at the surface of nucleus, in the 
potential Vq. Its kinetic energy is equal to Tp + Vq. Relativistic kinematics 
is used for description of the reaction (i.e. the total energy of a nucleon is 
connected with its momentum and mass by relation: E'^ = + m?). All 
nucleons are propagated in time; their momenta and positions are evolved in 
time as follows: r{t + 5t) = r{t) + p{t + 5t) = p{t) - VrVo6t. 
At time t = 0, the incident nucleon is hitting the nuclear surface. Next, all 
particles are moving along straight line trajectories, until two of them reach 
their minimum relative distance, or until one of them hits the nuclear surface. 
When a particle hits the nucleus surface from inside, two cases are consid- 
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ered. If the kinetic energy of the particle is lower than the Vq, it is reflected 
on the surface. If the kinetic energy is higher than the Vq, the particle is 
transmitted randomly with some probability, see Ref. [T7] . 
A collision takes place, if the minimum relative distance dmin between two 



particles fulfills the condition: dmin < y "^^"'^^^-^ , where: ^/s is the energy in 
the center of mass of the two particles, atot is the total collision cross section. 
The two particles are scattered elastically or inelastically, according to the 
energy - momentum conservation. Inelastic collisions with high probability 
lead to the formation of Deltas (A's). Mass of A is introduced with the 
Lorentzian distribution centered on the mean value equal to 1232 MeV, with 
the width 1=110 MeV. 

The following reactions are considered in the model: 
A^A^ ^ NN, NN ^ NA, NA ^ NA, AA ^ AA, A ^ Nn. 
All used cross sections and angular distributions are based on available ex- 
perimental data [17] . 

Additional condition restricting a collision is Pauli blocking. Collisions and 
A - decays are avoided, when the presumed final states are already occupied. 
The final states phase - space occupation probability (/j and fj, where i and 
j denote two particles predicted to be created in the final state) is evaluated 
by counting particles of the same kind, inside a reference volume in phase - 
space. The collision or decay is realized, when: Pij = (1 — /j)(l — fj) is larger 
than a random number chosen between and 1. 

The interaction process is stopped at time t = tgtop, determined by the aver- 
age behaviour of some quantities (e.g. an excitation energy of the nucleus, 
see Ref. [17]). At the end of the cascade, all remaining A's are forced to 
decay. 

2.2 Quantum Molecular Dynamics model 

Quantum Molecular Dynamics (QMD) fl9[ [2l] model has been developed 
mainly in order to investigate fragments formation during proton - nucleus 
or nucleus - nucleus collisions. The model, as a N-body theory, describes the 
time evolution of correlations between particles, what is essential in consid- 
eration of the fragments formation. 

In the approach, nucleons are spread out in phase - space with a Gaussian 
distribution. The coordinates and momenta of nucleons are designated si- 
multaneously. The wave function of each nucleon (i) is assumed to have the 
following Gaussian form: 





(2.1) 
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2.2 Quantum Molecular Dynamics model 



where: the Gaussian width - L = 1.08/m^ - corresponds to a root mean 
square radius of the nucleon of 1.8 fm, fio are the coordinates of the centers 
of the Gaussian wave packets. 

The normahzed Gaussian function represents one nucleon. 
In the model, Szilard-Wigner densities are applied. Those give a semi- 
classical approximation, depending simultaneously on coordinates and mo- 
menta. The Szilard-Wigner density is defined by the following expression, 
constructed of a wave function: ip{xi, ...,Xn) [25] : 

P(xi, = (— )" / ... / dyi...dyni>{xi + yi, Xn + Vn)* 

ilTT J — oo J — CO 

^{xi - 2/1, x„ - 2/„)e^^^^^^+-+^"^"i2.2) 

This is called the probability function of the simultaneous values of coor- 
dinates: and momenta: pi,...,pn. It has the following properties. 
Integrated with respect to the p, it gives probabihties for the different values 
of the coordinates: 

P{Xi,...,Xn;Pl,...,Pn)dp= \i){Xi, ...,Xn)\'^. (2.3) 

Also, integrated with respect to the x, it gives quantum mechanical proba- 
bilities for the momenta pi, ...,p„: 

/oo /"OO _^ 

... / ?/'(xi,...,x„)e-(pi^'i+-+P"^")rfa;i...dx, 
-oo J—oo 

(2.4) 

Based on upper definition, the Szilard-Wigner representation of considered 
in the model system is given by: 

The boundary distributions, i.e. the densities in coordinate and momentum 
space are given by: 

n{R,t) = J f{R,p,t)d'p = E '2!°^^^^' ^' ^^-^^ 

g{p,t) = J f{R,p,t)d^R. (2.7) 

In order to construct the initial system, the centers of the Gaussians (i.e. 
nucleons) are chosen randomly in coordinate and momentum space, in the 
following way. First, the positions {fl) of nucleons are determined in a sphere 
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of the radius R = 1.14 ■ A^^^. The numbers are chosen randomly, rejecting 
those which would position the centers of two nucleons closer than r^m = 1-5 
fm. In the next step, the local density (p(?^)), at the centers of all nucleons, 
generated by all the other nucleons is determined. Then, the local Fermi 
momentum (pp) is calculated: pp = h(^7r'^p(fi)Y^^. Finally, the momenta 
(pi) of all nucleons are chosen randomly, between zero and the local Fermi 
momentum (p^). Then, all random numbers, which position two nucleons 
in phase - space closer than: {fl — fj^ipi — PjY = dmin are rejected, and 
must be chosen again. The initialization process lasts long time. Typically, 
only 1 of 50 000 initializations is accepted under the criteria. But the finally 
accepted configurations are quite stable, usually no nucleon escapes from the 
nucleus in 300 fm/c, see [TH]. 

After successful initialization, the nuclei (in case of nucleus - nucleus 
collisions) are boosted towards each other with the proper center of mass 
velocity. 

During propagation, only positions (r^) and momenta (pi) of nucleons (i) 
are changed, the width of the wave function is kept fixed. The mean values 
{rio,Pio) are evaluated at time under the influence of two- and three- body in- 
teractions (important for preserving the correlations and fluctuations among 
nucleons), according to the classical Newtonian equation of motion: 

_ dH 

dH , , 

P. ^ (2.8) 

with the Hamiltonian H = T + U, where: T is total kinetic energy and U is 
total potential energy of all nucleons. 

Usually, the differential equations are solved using an Eulerian integration 
routine with a fixed time step At, where the momentum is evaluated at time 
points halfway between the times of the position determinations: 

Pio{n + 1) = Pio{n) - VrUi{n + ^)At, 
r„(n + i) = r.„(„ - \) + ^-^l^L^^^t + V,C/,(„)A*. (2,9) 

Assumed in the model total interaction is composed of a short range 
interactions between nucleons {V^"'^ ~ 5{ri — r2) + 5{ri — r2)5{ri — rs)), a 
long range Yukawa interaction (y^'^^ ^ ^\rl-r2\/i' ' ^ ~ ^ charge 

Coulomb interaction {y'^°'^^ ^ in^rl] ' ^1-2 stands for charges). 
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Potential acting on each particle is given by the expectation value of the two 
and three body interaction: 

U,it) = uPit) + uPit) (2.10) 

(2) 

where the two body potential is given by local, Yukawa and Coulomb 

(3) 

interaction terms, while the three body potential t/j includes only the local 
interactions term, see Ref. pTU] . 

Two nucleons can collide if they come closer than r = yja/n, where a is 
a total nucleon - nucleon cross section. Additionally, Pauli principle is taken 
into account. In the model, the measured free nucleon - nucleon scattering 
cross section is used. However, the effective cross section is smaller because 
of the Pauli blocking of the final state. It means, whenever a collision has 
occurred, the phase - space around the final states of the scattering part- 
ners is checked. It is calculated, which percentage. Pi and P2, of the final 
phase - space for each of the two scattering partners, respectively, is already 
occupied by other nucleons. Then, the collision is blocked with a probabil- 
ity Puock = P1P2, or allowed with the probability 1 — Puock- If a collision 
is blocked, the momenta of scattering partners are kept with values, which 
they had before scattering. 

The scattering angles of the single nucleon - nucleon collisions are chosen 
randomly, with the probability distribution known from the experimental 
nucleon - nucleon scattering |26j . 

Inelastic collisions lead to the formation of Deltas, which can be reabsorbed 
by the inverse reactions. 

It is assumed in the QMD approach, that only these beam energies are ac- 
cepted, at which no more than 84% of all collisions are blocked. Therefore, 
the low energy limit of the model is kinetic energy Tiab = 20 MeV/nucleon 
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2.3 Percolation model 

Percolation model has been introduced in the eighties as minimum informa- 
tion approach, based on purely topological and statistical concepts. It was 
in order to describe fragment size distributions, as an outcome of nuclear 
fragmentation process within the simplest possible physical framework |27] . 
However the model is flexible enough to allow for the inclusion of different 
physical mechanisms. 

In general, one distinguishes between bond and site percolation model. 

Site percolation corresponds to the case where each site is either occupied 
by one particle or is empty. One fixes a probability p e [0, 1], and for each 
site generates a random number ^, taken from a uniform distribution in the 
interval [0, 1]. If ^ < p the site is said to be occupied, if ^ > p it is empty. 
Checking all sites, one obtains an ensemble of occupied and empty sites. Sets 
of occupied sites are called clusters. Finally, space is topologically covered 
with clusters and empty space. 

Bond percolation corresponds to the case where each site is occupied, 
but neighbouring sites being bound or not to each other. The number of 
neighbouring sites is fixed by the bonds, which depend on the geometric 
structure of space occupation. A similar procedure, as in case of site perco- 
lation, works also in this case. For a fixed bond probability p one considers 
a pair of neighbouring sites and generates a random number ^ G [0,1]. If 
^ < p the sites are linked by a bond, if ^ > p they are not. Checking all 
possible bonds between neighbouring sites one obtains clusters made of con- 
nected particles. There are clusters of different sizes, appearing with a given 
multiplicity, which depends on p. 

The first applications of percolation concepts were proposed by W. Bauer 
and collaborators [2HI EH] , Campi and Desbois [30] . While Campi and Des- 
bois used a site percolation, Bauer with his group have proposed a model 
based on bond percolation theory. Nowadays, the bond percolation has been 
used as a base for models of fragmentation more often than the site one. 
Below, an outline of model proposed by Bauer et al. is presented. 
As any percolation model, it is based on two crucial ingredients: a descrip- 
tion of the distribution of a set of points (i.e. nucleons) in a space and a 
criterion for deciding whether two given points are connected. The target 
nucleons are represented by points occupying uniformly an approximately 
spherical volume on a simple cubic three-dimensional lattice in coordinate 
space. The lattice spacing d is computed from the normal nuclear density: 
d = ~ 1.8 fm, where: po = 0.16 nucleons/ fm^. The number of points is 



equal to the number of target nucleons and is conserved during the calcula- 
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tion, therefore the conservation law of mass in the calculated fragmentation 
process is fulfilled. Initially, each nucleon is connected by bonds (representing 
the short-ranged nuclear interactions) to its maximum six nearest neighbors, 
depending on its location in the target. Then, a point-like proton, with an 
interaction radius r colhdes the target. Because it is assumed that the mo- 
tion of the nucleons in the target is neglected, the projectile sees a frozen 
image of the target (it is feasible, as typical Fermi motion speed of nucleon 
is significantly lower than speed of incoming proton). For a given impact 
parameter b, the proton removes from the lattice nucleons occupying a cylin- 
drical channel with the radius r along his straight path in the target. It 
is typically 6-8 nucleons. All of remained nucleons, called spectators, are 
still connected via bonds. These bonds are then broken with a probability 
p, which is a percolation parameter. The parameter p should be related to 
some physical input. For example, it is reasonable to assume that p is a linear 
function of kinetic energy of the projectile, and also an increasing function 
of the excitation energy (E*) of the spectators: p = E*/Eb, where: Eb is 
the nuclear matter binding energy per nucleon (16 MeV) [HT] . 
The breaking probability has to be also dependent on the impact parame- 
ter b of the proton. Bauer et al. have used following dependence: p{b) = 
i+exp[{b-R)/a] ' ^^cre: Po = p{b = 0), i? is a radius of the target nucleus, 
a = 1.0 fm is a diffuseness parameter |29l[3T]. 

For given parameter b, the breaking probability is assumed to be uniform for 
all bonds, independent of their position on the lattice. 

Using the breaking probability p as an input parameter, a Monte-Carlo al- 
gorithm decides for each bond individually whether it is broken or not, as 
follows. For a given p{b), a random number C,ijk ^ [0, 1] is generated for each 
bond Bijk, where: the indices ijk correspond to the spatial location of the 
center of the bond on the lattice. If ^ijk > p, the bond Bij^. is unbroken, if 
^ijk < P, the bond is broken. Then a cluster search algorithm [22] is used 
to find out which nucleons are still connected by bonds i.e. form clusters. 
Taking into account all impact parameters, inclusive mass and multiplicity 
distributions can be obtained. That can be compared to experimental re- 
sults. 

It is surprising that using only one free parameter and simple geometrical con- 
siderations, this model is able to reproduce experimental mass yield curves 
with a good accuracy, in particular, the power law behavior (~ Ap^) at small 
masses and the U-shape distribution of the whole mass range. 

Bauer's group have used such model also to study the possibihty of observ- 
ing a phase transition of nuclear matter in collisions of high energy protons 
(> 10 GeV) with heavy targets. Since inclusive fragment mass distributions 
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follow a power law behavior: a{Ap) ^ Ap , for Ap < At/3, (A = 2.6 is 
independent of target mass for heavy targets) [32] similar to the mass yield 
distribution of droplets condensing at the critical point in a van der Waals 
gas ( (T(m) ^ m""^, with the critical exponent r = | ), they suggested that 
nuclear multifragmentation proceeds via a liquid-gas phase transition of nu- 
clear matter. 

Bauer et al. have accented the importance of their result for experimental 
study of the phase transition, that the critical events are not the ones with 
the highest multiplicities, but the ones with the highest value of standard 
deviation of mass distribution. 

In the percolation models, such as described above, it is assumed that 
nucleons are distributed uniformly in the sphere, and the total excitation 
energy is assumed to be uniformly distributed over the whole excited sys- 
tem of the spectators, what is equivalent to assumption that equilibrium is 
reached. In this picture, the angular distribution of fragments should be for- 
ward peaked, because the momentum transfer from incident proton to the 
nucleus is in average in forward direction. In reality, as the incident energy 
increases, the mass fragments angular distribution grows from forward to 
sideward {Ep > 10 GeV) [33] or backward {Ep > 100 GeV) j31] peaked in 
laboratory frame, what contradicts to the picture of the fragmentation from 
an equilibrated system, and indicates that the nucleons may not be uniformly 
distributed spatially and the excitation may depend on the position inside 
the nucleus. 

For the understanding of this sideward emission (for which not satisfac- 
tory explanation has been given so far) Hirata et al. [35] investigated the 
non-equilibrium dynamical effects, such as non-spherical nuclear formation. 
They formulated Non-Equilibrium Percolation (NEP) model, which they use 
in a combined framework with a transport model [33]. The main differences 
between equilibrium percolation and the NEP model are following. The ini- 
tial conditions of percolation, instead of putting nucleons on sites with some 
assumed occupation probability, is taken from the results of the dynamical 
transport model calculations. The bond breaking probability is assumed to 
be dependent on the position and momenta of the nucleons. It is calculated 
by considering excitation energy, distance and momentum difference between 
the nucleon pairs, instead of giving a common breaking probability for the 
bonds connecting nearest neighbor sites. 

Analysing fragmentation process with the NEP model, considering cal- 
culated effects, Hirata et al. have found following mechanism of sideward 
enhanced fragments emission. Based on fragments formation point distri- 
butions, both in case of central and peripheral collisions, the fragments are 
formed mainly near the surface of the nuclei. It is due to fact, that along the 
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incident proton path, nucleons collide with the leading proton or secondary 
cascade particles. Since they have large kinetic energies, they increase the 
bond breaking probability. As a result the fragment formation is suppressed 
along the incident proton path. Fragments are formed mainly in the cold 
region around the hot zone. Their formation points are distributed non- 
spherically, in a doughnut shaped region. Hirata et al. noticed that this 
effect alone does not generate any anisotropy in angular distributions. 
Based on analysis of fragments energy distributions, they found that the en- 
ergy distribution is effected by Coulomb repulsion. Calculated results after 
the Coulomb expansion well reproduce the qualitative behavior of experi- 
mental data [35j. This is Coulomb repulsion between formed fragments that 
pushes and accelerates them sideways of the doughnut region. It means, the 
Coulomb repulsion modifyies the angular distribution from forward peaked 
to sideward peaked. 

Bond percolation model with non-equilibrium effects, have been investi- 
gated also by Yamaguchi and Ohnishi |36j. They introduced to the model 
isospin dependence, by assuming that neutron-neutron and proton-proton 
bonds are always broken, while neutron-proton bonds make the nucleus 
bound. Additionally, by comparing calculated fragments energy spectra with 
experimental data, they have got an agreement, when taken density of frag- 
menting nuclei is not equal the normal nuclear density po, but of around 
~ Po/3. It means, target nucleus, after being heated by incident proton, 
would expand up to rather low, mechanically unstable density. 
They have also considered, that the incident proton heats up either cylindri- 
cal or conic shaped region around its path in the target. They have found, 
that in case of cylindrical heated region, formed fragments are pushed, by 
Coulomb repulsion, more strongly in sideward directions. If a conic shaped 
region is heated up, fragments are pushed in rather backward directions. 
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Specific models for fast stage of 
proton-nucleus collision 

In the frame of this work, dynamical analysis of fast stage of proton - nucleus 
reactions are performed within transport approaches: Boltzmann-Uehling- 
Uhlenbeck (BUU) [IHIEZ! and Hadron String Dynamics (HSD) [SHI El] mod- 
els. The models have been specially developed in order to enable description 
of considered here reactions. 

The BUU model is used to calculate proton - nucleus reactions in projec- 
tile kinetic energy range only up to about 2.5 GeV. That is because included 
in the model processes go mainly through single resonances excitations (i.e. 
A, A^(1440), A^(1535)), what is correct in this energy range. While an in- 
cident energy increases, the density of produced resonances also increases. 
In this possible proper description of processes requires taking into 

consideration hadron - hadron reactions on the level of elementary quark 
- quark interactions. This can be done by employment of a string model 
(e.g. FRITIOF model [IQ]), where during inelastic collision, two interact- 
ing hadrons are excited due to longitudinal energy- momentum transfer [lO] . 
The formed excitation, so-called string, represents a prehadronic stage. It is 
characterized by the incoming quarks and a tubelike colour force field [H] 
spanned in between. The string is then allowed to decay into final state 
hadrons, with conservation of the four-momentum, according to e.g. Lund 
string fragmentation model ^2] . 

The HSD approach includes the FRITIOF scheme of string dynamics and 
the Lund model of hadron production through string fragmentation. It is 
employed, in particular, for incident energies higher than projectile energy 
2.5 GeV. For projectile energies lower than about 2.5 GeV, the BUU code is 
used in the frame of the HSD model. 
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3.1 Transport equation 



Below, a description of the approaches is given. 

Both the models are based on transport equation. 

3.1 Transport equation 

Historically, the transport equation originate from classical Boltzmann equa- 
tion for one-body phase-space distribution function f{f,v,t) defined such 
that 

/ /(f, V, t) d^v (Pr is the number of particles at time t positioned in element 
volume (Pr around r, which have velocities in volume element of velocity 

space d^v around v. 

Let's consider particles, in which an external force F with mass m acts and 
assume initially that no collisions take place between the particles. In time 
t + 5t the velocity v of each particle will change to v + {^)5t and its position 
r will change to r -|- vSt. Thus the number of particles f{f,v,t) d^vd^r is 
equal to the number of particles f{f+ vSt , v + {—)St , t + dt) d^r'd^v', what 
is explained by the Liouville theorem: 

The volume of phase-space element is constant, if movement of 
all particles inside is consistent with canonical Hamilton equation 
of motion. 

and written as: 

F 

f{f+ vSt,v + {—)St ,t + St)- f{f, v,t)^0 (3.1) 

If collisions occur between the particles, an additional element, i.e. colli- 
sion term is needed. This gives the following equation describing evolution 
of the distribution function: 

/(f + vSt,v+ {^)St ,t + St)- f{f, V, t) = {df/dtUiSt (3.2) 

Letting St ^ and expanding into the Taylor series gives the Boltzmann 
equation: 

— * 

iid/dt)+Wr ■ v + {-)Vv)fir,v,t) = (df/dtUi (3.3) 

An apparent form of the collision term {df /dt)coii can be found considering 
an element volume A at time t, around position (r , v) and an element volume 
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B at time t + 5t, around position (r + v6t,v + {^)6t). These two element 
volumes are so similar, that letting 5t — > 0, particles knocked out from A, 
due to coUisions, will not get into B. Particles being outside A, during time 
5t, will get into A, and they will be inside B. So, the number of particles 
inside B, at time t + 6t, at dt ^ 0, is equal to the initial number of particles 
inside A, at time t, and a magnitude of relative modification of number of 
particles due to coUisions, during time 5t. Therefore, a form of the coUision 
term {df/dt)coii can be calculated as a difference between the number of 
collision in a time range {t,t + St), when one of particles after collision is 
situated in element volume d^rd^v around position {f,v), and the number of 
collision in a time range {t,t + St), when one of particles before collision is 
situated in the same element volume d^rd^v around position (r , v) . It can 
be done by assuming that the density of particles is low enough, that only 
binary collisions need be considered. It is also assumed that the velocity 
of particle is uncorrelated with its position in the space. It means that in 
element volume d^r the number of particles pairs with velocities in volume 
elements of velocity space d^Vi around Vi and d^V2 around V2 is equal to: 
[/(f, v-{, t) d^r d^V,] [/(f, V2, t) d\ d^V2] . 

The number of binary collisions {vi, V2 V3, vl) inside element d^r, in 
time range {t , t + St) is equal to: [/(r, V2, t)d^V2]\v[ — V2\a{fl)dil,St, 
where: 

vi and V2 are the velocities of the two particles before collision, 
?J3 and vl arc their velocities after the collision, 

a{Q) is the differential cross section for a reaction, in the centre of mass ref- 
erence frame, 

is the sohd angle the particles are scattered into (the angle between vectors 

\vi — V2 \ is the magnitude of the particles relative velocity before the collision, 
[f{r,V2,t)d^V2]\vi — 1^2! is the density of particles flux equal to the product 
of particles density and their velocity. 

The total number of coUisions, where one of the particle before coUision is 
situated inside element d^rd^Vi around {r,v[) is obtained multiplying the 
number of binary collisions by number of particles with velocity vi, inside 
element d^r and integrating over all possible V2 and O: 

j d^V2 J dna{fl)f{f, V2, t)\v{ - V2\[f{f, v{,t)dh d''vi]St (3.4) 

Taking into consideration the inverse binary collision: {ih,, ^ Vi, 7)2), and 
using analogical method as above, the total number of collisions, where one 
of the particle after collision is situated inside element d^rd^vi around (r, to{) 
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3.1 Transport equation 



is obtain: 

'rfV / dQa\Q)f{r,v7,t)\v3-v,\[f{f,vs,t)d\d\3]St (3.5) 



Because collisions vl, V2 ^ vs, vl and vs, vl ^ vl, V2 are inverse collisions, 
so: a{n) = a'{n). 

From covservarion law for energy and momentum: jt'i — i'2| = |'^~'^|- 
From Liouville theorem: d^vid^V2 = d^v^d^v^. Subtracting equations (I3.5p 
and (13.41) . and using above assumptions the collision term can be written as: 



{df/dtUi = J d'v2d'vsa{n)dn\vi - V2\{hh - /1/2) (3.6) 

where: /i = /(f, ^Ti, t), /a = f{f,V2,t), fa = f{r,V3,t),f4 = f{f,V4,t). 

Joining equations (13. 3p and (13.61) one obtains the classical Boltzmann 
equation: 

+ ( A + dU{r,put) ^ d_ _ dU{f,put) t) = 

dt mi dpi df df dpi ' ' 

J d^p2d'p3dQ a{Q)\v-l - V2\ 6\pi+p2-P3-P4) ■ [fsfA - /l/2](3.7) 



{2nY 

where: 

V = p/m, 

F = —{dU{r)/df), Uif) is position dependent potential. 

In 1933 Uehling and Uhlenbeck have developed the equation by adding 
the Pauli factors ^3] . 

Due to Pauli blocking, a collision can occur only if in the final state, there 
are free quantum states. 

Probability of finding the free quantum state in a phase-space volume is equal 

to: [1 — /(r, p, t)], what is responsible for fermion Pauli blocking. 

The probability of two particles collision with momenta pi and p2 is equal 

to: 

[1 - /(r,P3,t)][l - /(r,P4,t)] = 73/4 (3.8) 
Probability of inverse occurrence is equal to: 

[1 - /(r,Pi,t)][l - /(r,P2,t)] = IJ2 (3.9) 
Expression (13. 8p and (13.90 are Pauli factors. 
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Including in fl3.7p the Pauli factors gives: 

+ ( A + 9U{r,Put) d_ _ dU{f,pi,t) t) = 

dt uii dpi dr dr dpi 

4 f 

I d^p2 d^P3 dQ a{n) \vi - V2\ S^{Pi + P2 - P3 - Pi) 

■[hUlil2-hf2lJ,] (3.10) 



(2vr) 

The equation is named Boltzmann-Uehhng-Uhlenbeck (BUU) equation. 



3.2 Boltzmann-Uehling-Uhlenbeck model 

The theory based on the transport equation (13.101) (it means the Boltzmann 
equation with a self-consistent potential field, and with a collision term that 
respects the Pauli principle) was used first time to nuclear collisions descrip- 
tion by Bertsch, in 1984 [44]. 

The BUU equation is solved numericaly, using Monte Carlo method, repre- 
senting the one-body phase-space distribution by discretized test particles: 

. N-A{t) 

fir,p,t) = - Y: 5\t - nit)) 5\p - m) (3.11) 

where: 

is a number of test particles, 
A{t) is a number of real particles at time t 



Likewise, one collision is replaced by parallel collisions. 



All of the test particles give part to the density of nuclear matter not in 
single points, but they are smeared with Gauss distribution. This way, the 
effect of quantum smearing is included. 
The density is calculated on the grid r^: 



1^ 1 .-fr;-n^2 



where: A is the Gauss width parameter (taken usually equal 1). 



The initial coordinates of particles of target nucleus have Wood-Saxon 
distribution form: 

Po 

1 + exp{^) 



P(0 = ^-^^ (3.13) 
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where: 

a = 0.025 ■ ^ 0.29 [fm] = 0.5 [fm] 

R = 1.124 ■ A^/^ [j^] 

Po = 0.168 [nukleon/ fm^]. 



Initially, the target nucleus is in the rest, the total momentum (i.e. the 
sum of momenta of all particles) is equal to 0, but the local momenta of 
particles are determined homogeneously on the Fermi sphere with radius 
PF{r): 

PFir) = (^p(r))i/3 (3.14) 

The projectile is a single proton. The test particles replacing the proton 
are distributed homogeneously on a thin cylinder, with a radius equal the 
radius of target. Thanks to this approach, each test particle has different 
impact parameter, so the results of calculations are averaged over all impact 
parameters. 

The solution of the transport equation is the single-particle phase-space 
distribution function, depending on time. 
The collision is numericaly evolved by fixed time steps. 
The test particles propagate between collisions according to the classical 
Hamilton equations of motion: 

A dU{fi,pi,t) 

P. = 9^ (3-15) 

dU{fi,pi,t) 



Ti = pi/ym"^ + p^ -\ — (3.16) 

where: U{ri,pi, t) is a mean field potential, dynamically changing, calculated 
function of local density: 



U{r) = hop{f) + -Up{^)'/' + Vof dV'^-^^P^-^^^^pi?) + Vcoui (3.17) 



where: to=-1124 MeV-fm^, ^3=2037 MeV■fm^ ^=-378 MeV, p=2.175 fm'^, 
see Ref. 



The momentum and coordinates of all particles taking part in reaction 
are calculated in the successive time steps. Using the values of momentum 
and coordinates, all other quantities (i.e. nucleon density, mean field poten- 
tial) are calculated. 
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The nucleoli - nucleoli collision at a fixed time step is introduced as fol- 
lows: when two nucleons come closer than the distance bmax = \[^WifJ^ 
(where cr]^/lf is the maximal cross section for nucleon - nucleon interaction in 
nuclear matter (30 mb [H])) they are made to scatter, but if the final state 
is Pauli blocked, this collision is canceled. 

The BUU model describes the propagation and mutual interaction of nu- 
cleons, Delta's, A^* - resonances, and also vr and t] - mesons. 
In the model, the following reaction channels are included: 

A^A^ ^ NN, 
NN ^ NR, 
NR ^ NR, 
AA ^ NR, 
R ^ Ntt, 
A^(1535) ^ Nr], 
NN ^ NNtt, 
ttN 7tN, 
ttN —>■ ttttN, 

where: is a nucleon, R stands for a resonance A, A^(1440) or A^(1535). 
Cross sections for the reactions, used in the model calculations, are parametriza- 
tions of the experimental cross sections taken from [45]. 

All resonances are allowed to decay into two particles (besides the decays 
due to collisions with other particles, e.g. NR NN). The decay of a 
resonance is determined by its width r(M). The decay probability P is 
calculated in every step of time, according to the exponential decay law: 
P = 1 — e~^^^'^'>^^^"''>^^ , where: T{M) is the energy dependent width of the 
resonance, 7 is a Lorentz factor related to the velocity of the resonance, 6t is 
a time step size of the calculations. 

For the A decay the parametrization given by Koch et al. [16] is used. Details 
concerning higher resonances can be found in [57]. 

In each step of time it is decided if the resonance may decay and to which 
final state it may go. If the final state is Pauli blocked, the resonance decay 
is rejected. 

3.3 Hadron String Dynamics model 

The HSD model is based on the same transport equation fl3.1UI) as the BUU 

approach and solved also by use of the test particle method. 

In the HSD model, propagation of the following real particles are included: 
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baryons (p, n, A, A^(1440), A^(1535), A, S, S*, H, Q), the corresponding 
antibaryons and mesons (tt, i^T, ?], r]', p, u, 0, if*, ai). 



3.3.1 Cross sections 

The low-energy baryon-baryon and meson-baryon colhsions (i.e. with 
the invariant energies below "string threshold" : ^/s < 2.65 GeV and ^/s < 2.1 
GeV [Sn], respectively) are described using the explicit cross section, as in 
the BUU code. The following parametrization of the experimental total and 
elastic p + p, p + n, vr+p, 7i~p, K~^p, K~p cross sections, taken from |15] , is 
used: 

a{p) = A + B-p"" + C-ln^{p) + D-ln{p)[mh], (3.18) 

where: p is a momentum of incident proton in laboratory frame, 

A, _B, n, C, D are constants, with adequate values for different processes 

cross sections, see Ref. ^8] . 

For reaction channels: ttN — > vrA^ and ttN — >■ ttttN, to have them consistent 
with the experimental inelastic pion-proton cross section below the string 
threshold, instead of fl3.18p . the following parametrization is used |l8] : 

atotiPn) = 25.0 + lQ.O-p~'-^[mb], 

aeiip^) = 3.5 + 14.0 ■ p~' '[mb], (3.19) 

where: is a momentum of incident pion (in GeV/c). 
Additionally, the channels: 
pN Ntttt, pA —>■ Ntttc and 
ujN —y Ntttttt, ujA Ntttttt, 

are included with an energy independent cross section equals to 30 mb. 
Because of very low number of produced hyperons in low-energy proton- 
nucleus collisions, the hyperon(A, S) - nucleon interactions are neglected in 
the model. 



Angular distribution of elastic collisions depends on energy |18] , therefore 
the following parametrization of the differential elastic nucleon - nucleon 
cross section, taken from Cugnon et al. [26J, is applied: 

where: s is the invariant energy of collision squared (in GeV), 
t is the four-momentum transfer squared, 

X _ r. [3.65(v/i-1.8776)]'^ 
^y^) — "l+[3.65(v^-1.8776)]6- 
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In the approach, the high-energy elastic and total baryon - baryon 
and meson - baryon collisions (i.e. with energies above the "string thresh- 
old"), are related to the measured cross sections by: 

^f{V~s) = ofi^s) = atn^s) = ... = 0.5«;^(Vi) + <r^(Vi)), 

^r(v^) = 4^(v^) = - = <^T'iv~s), 
^er(x/i) = <^^i^iv~s) = ... = 

where dots stand for other combination in the incoming channel, with K — 
{K+, K°, K*+, K*^) and K = {K-,K^, K*' ,K*^) . The same relations are 
applied for the total cross sections. 

The high-energy inelastic baryon - baryon and meson - baryon cross 
sections obtained by this procedure are equal to 30 mb and 20 mb, re- 
spectively. This corresponds to the typical geometrical cross section (i.e. 
'^inei ~ Tri?^) , SO it should be reasonable input for the calculations. 
Due to this fact, in order to include all baryon - baryon and meson - baryon 
high-energy inelastic cross sections, only final state rates must be specified. 
Therefore, the string model (i.e. FRITIOF model) is employed, which will 
be described hereafter in this work. 



In the HSD model, the following meson - meson reaction channels are 
included: 

7rp <-> 0, 
7rp <-> ai, 

KK ^ 0, 
TiK ^ K*. 

They are more probable at high-energy proton - nucleus collisions, but as 
they appear as secondary or higher order reactions, the average energy for 
such processes is rather low. In this case, the cross section within the Breit- 
Wigner parametrization is employed. Therefore, the reactions a + b ^ rriR 
— J. c + (i, where: a, 6, c, d are the mesons in the initial and final state, 
respectively, and uir denotes the intermediate mesonic resonance (p, Oi, 0, 
K*), are described by: 

where: 

Sa and Sb are spins of the particles, Jr is spin of resonance; 

Tn^ab and T^^cd are partial decay widths in the initial and final channels, 

Mr is the mass of the resonance, Ttot is the total resonance width. 
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Pi is the initial momentum in the resonance rest frame. 
The decay widths and the branching ratios for the mesonic channels are 
adopted from the nuclear data tables [15], without introducing new param- 
eters. 

Additionally, strangeness production in meson - meson collisions is included, 
with an isospin averaged cross section |49|: 




where: Sq = 4:mj^, mm stands for all possible non-strange mesons in the 
incoming channel (e.g. nir —>■ KK, Tip — > KK). 



Quarks, as colour- charged particles cannot be found individually. They are 
confined in two- or three-quarks systems, i.e. colour neutral hadrons. The 
quarks, in a given hadron, exchange gluons. If one of the quarks is pulled 
away from the other quarks in the hadron, the colour force field, which con- 
sists of gluons holding the quarks together, stretches between this quark and 
its neighbers. Because of interaction between gluons [41j, the colour field 
lines are not spread out over all space, as the electromagnetic field lines 
do, but they are constrained to a thin tube-like region. While the quarks 
are pulled apart, more and more energy is added to the colour force field. 
So, in such a formation, the four-momentum can be accumulated. At some 
point, it is energetically possible for the field to break into new quarks. The 
four-momentum is conserved, because the energy of the colour force field is 
converted into the mass of the new quarks. Finally, the colour force field 
comes back to an unstretched state. 

With this picture in mind, high-energy hadron - hadron interactions models, 
called string models [39] have been created, where the formation composed 
of quarks and the colour field between is called a string. 

In the HSD approach, in order to describe the high-energy inelastic hadron 
- hadron collisions, FRITIOF model is applied [ID]. In the model, hadronic 
collision corresponds to large longitudinal and small transversal energy - mo- 
mentum transfer. It means, to the stretching of longitudinally extended 
string-like colour force field along the beam direction, between constituent 
quarks of the incoming hadron. The created excitation, i.e. string is a dy- 
namical object, which may decay into final state hadrons, according to the 
Lund fragmentation scheme [32], implemented in the FRITIOF model. The 
field between quarks is confined into a tube, called "flux tube", which is 





(3.22) 



3.3.2 String Model 
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one-dimensional object. The uniform colour field contains constant amount 
of energy stored per unit length. The total energy {E) of the field is pro- 
portional to the length (L): E = n ■ L, where k = 1 GeV/fm, is a string 
tension ^2]. Due to longitudinal energy- momentum transfer, spread over 
some region, the colour separation occurs, i.e. as seen in CM string frame, 
there will be two extended parts of the string moving forward and backward, 
along the beam direction. The potential between quarks is linearly rising. 
As result, the system breaks. Because of colour confinement there is never 
a single quark in isolation. After the string breaks, on its ends new quarks 
appear. The new qq pair is created from the available field energy. By that 
means, the energy of the initial string decrease (a part of the energy is used 
for qq pair production), but the breaking process is not finished yet. The 
quarks of the new created strings are also moving in opposite directions in 
the strings rest frames. Thereby the original system breaks into smaller and 
smaller pieces, until only physical hadrons remain, (i.e. baryons as bound 
systems of three quarks, antibaryons - three antiquarks and mesons as quark 
- antiquark systems). 

In the HSD approach, baryonic {qq — q) and mesonic {q — q) strings are con- 
sidered. 

The implementation of the string model into the transport approach implies 
introduction of a time scale for the particle production processes. The time 
scale is given by a formation time tf, which includes formation of a string, 
fragmentation of the string into small substrings due to qq and qqqq pro- 
duction, and formation of physical hadrons. The formation time should be 
also related to the spatial extension of the interacting hadrons. It means, it 
should be big enough so that quark - antiquark pair could reach a distance 
corresponding to a typical hadron radius, i.e. 0.6 - 0.8 fm |3H|- In the HSD 
model, the formation time is a single fixed parameter for all hadrons, it is 
set to tf = 0.8 fm/c in the rest frame of the new produced particles |50j . 
The particles production proceeds as follows. 

If a system contains originally e.g. q^ and qo, moving in opposite directions 
with large energies, it breaks after some time into two parts by a production 
of a pair qi g^^, at a space-time point (xi, ti). The new produced quarks also 
move in opposite directions. Two subsystems are created by TJq qi and q^ go- 
As a result, a colour force field between the new pair qi q^ vanishes. At a 
later time another pair q2 ^2 can be produced at (x2, ^2), due to breaking e.g. 
qi go subsystem. Analogically, new subsystems g^ g2 and g2 go are created. 
A colour force field between the new quark pair g2 g2 vanishes. Created sub- 
systems either are hadrons or else will fragment further, until only hadrons 
remain. In the model, string fragmentation starts always in its center (in the 
CM frame of the string). All quark pairs production points are separated 
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in a space-time. A total momentum in the rest frame of a string is equal 
to zero. The energy and momentum in the production process is conserved. 
Masses of the finally produced hadrons are equal to the masses of physical 
hadrons. 

In the HSD model, the production probability (P) of massive ss or qqqq 
pairs is supressed in comparison to light quarks pairs production {uu, dd). 
Inserting the following constituent quark masses: = = 0.3 GeV and 
m, = 0.5 GeV, one gets: = 0.3 and ^^Mi = o.07. 

* ' F(uu) F(uu) 

So, the suppression factors used in the model are: 
u : : s : = 1 : 1 : 0.3 : 0.07, see Ref. [IH]. 
As a result, mainly mesons are produced. 



The model assumes that there is no final state interaction of the produced 
hadrons included in the model. 

Because most of the strings, in a given space-time volume, fragment within 
a small time interval, the interaction of the string field spanned between the 
constituent quarks with other hadrons is not taken into account. But the 
secondary interactions of the quarks or diquarks inside the strings are con- 
sidered in the approach. The following cross sections of such interactions are 
used HH]: 

a{q -B) = \a{B - 5) ^ 10 mb, 
a\qq - B) = la{B - B) ^ 20 mb, 
a{qq — q) = — B) ^ 6.6 mb. 

Because most of strings are stretched longitudinally, parallelly to each other, 
and radius of string is small, equal to about 0.2 - 0.3 fm |5T], there is no 
string - string interaction included in the HSD model. 



The characterized above version of the HSD code has been used and 
developed in order to describe proton induced spallation reactions in wide 
energy range, and mass of target nuclei in the following aspects. First of 
all, to calculate properties of residual nuclei remaining after first stage of the 
reaction, forming an input for models of the second stage calculations. 
Additionally, the code has been developed for description of pions produced 
in proton - nucleus spallation reactions. Moreover, a version of HSD code 
that allows for calculations of pion induced reactions has been prepared. 
Each of the directions of the HSD code development will be presented and 
discussed in the next sections of this work. 
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3.4 Further development of the model 

As a result of the first stage of proton - nucleus reaction, apart from emitted 
particles (mainly nucleons and pions), a hot excited nucleus remained. Prop- 
erties of the residual nucleus, i.e. mass (Ar), charge (^r), excitation energy 
(E^), three - momentum (pji) and angular momentum (Lji) form an input 
for the second stage models. In order to calculate the properties, the original 
version of the HSD code has been modified. Corrections ensuring energy and 
charge conservations have been made, which are important here, but played 
minor role in the original version. Also, corrections enabling calculations of 
various quantities in function of time, for very large times of propagation 
have been found. 

The properties are evaluated in the following way. 

First, using that in the models calculations, four-momenta of all hadrons are 
propagated in time, the particles that have left the residual heavy fragment 
are identified. This can be done in two ways. In case of low-energy collisions, 
a sphere of observation, with radius equal to Ra + 2 fm, where: Ra denotes 
the radius of the target with mass number A, can be considered. All the 
particles inside the sphere are treated as belonging to residual nucleus, the 
particles outside - as emitted. In case of high-energy collisions, when target 
nucleus is moving faster during reaction, it is easier to consider a baryon 
density criterion. It means, particles, which are positioned in baryon density 
lower than 0.02 nucleon/fm^ are treated as emitted, the rest of the particles 
form a residual nucleus. At intermediate proton impact energy (0.1 - 2.0 
GeV) both methods are equivalent. In the frame of this work, the second 
criterion is used. Nevertheless, both of the methods of classification are con- 
nected with some inaccuracy, concerning the very low energetical particles. 
Distributions of particles escaped from the sphere of observation characterize 
absence of very low energetical part. In contrary, distributions of particles 
classified on the base of the baryon density criterion characterize overabun- 
dance in the very low energetical part. In this second case, it is because the 
density condition classifies incorrectly nucleons placed in the lowest density 
level of a residual nucleus, as emitted. Therefore, the density method needs 
to be completed by an additional condition concerning the particles kinetic 
energy. Particles in the nucleon density equal to 0.02 nucleon/fm'^ acquire 
momentum with a value from the range from zero to the local Fermi mo- 
mentum, equal to about 130 MeV/c, what corresponds to about 10 MeV of 
kinetic energy. Compliance of that information gives a complete and correct 
condition, i.e. particles in nuclear density lower than 0.02 nucleon/fm^ and 
with kinetic energy higher than 10 MeV are considered as emitted. 
Then, by exploring the conservation of total energy, mass number, momen- 
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turn and angular momentum, for each parallel ensemble (see description of 
models), the average values of properties of the residual nuclei are calculated 
as function of time, according to formulas: 

Npit) 

<E*j,>{t) = Etot- E ^pl + m^-Mn-Ec, 

i=l 

<AR>{t) = AT + Ap-Np{t), 

Np{t) 

<PR> {t) = Plot - E Pi(t)^ 
i=l 

<L>{t) = Ltot - E X Piit)^ (3-23) 

i=l 

where: Np(t) denotes the number of emitted particles, M/j is the mass of the 
residual nucleus. At is a mass of original target, Ap = 1 stands for incoming 
proton and Ec is the the energy of Coulomb interaction between the emitted 
particles and the residual nucleus. 

3.5 Stopping time criteria for the first stage 
model calculations 

In order to describe both stages of the reaction, model of the first stage must 
be used together with a statistical model for the second stage of reaction. 
The outcome of the transport model calculations determines an input for the 
second stage model. Thus, it is important to assume a proper duration time 
of the first stage calculations. In order to define the time at which the first 
stage calculations should be stopped, it must be verified, whether informa- 
tions obtained from the transport models are sensitive to the time duration of 
the first stage of the reaction. For this purpose, time variation of the average 
values of four physical quantities: excitation energy per nucleon, momentum 
in beam direction per nucleon, angular momentum and mass number of the 
excited residual nucleus after the first stage, have been analysed. As exam- 
ple, the dependences for p + Bi collision, at 3.0 GeV proton beam energy are 
discussed below. The average values of presented quantities are evaluated 
according to the equations (13.231) . 

First, let's look at the time evolution of the average value of excitation en- 
ergy per nucleon of the residual nucleus shown in Fig. 13.11 It is evident from 
the Figure, that it takes some time before proton approaches nucleus, i.e. 
about 10 fm/c (at the time zero, projectile is placed in some distance from 
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the target nucleus, ensuring that density distributions of projectile and tar- 
get are not overlapping). Next, up to about 18 fm/c, the excitation energy 
gains some maximal value, which corresponds to the energy introduced into 
target nucleus by incoming proton. Then, the excitation energy drops very 
quickly in the range 20 - 30 fm/c of the duration time of the first stage of 
the reaction, whereas it is varying only a little at larger times. It is seen that 
the average value of excitation energy per nucleon starts to stabilize at time 
35 - 40 fm/c. 

Looking at Fig. 13.21 it is seen, that the average value of the momentum of 




I ' ' ' ' 1 

10 20 30 40 50 60 

time (fm/c) 



Figure 3.1: Time variation of tlie average value of excitation energy per 
nucleon of the residual nucleus in p + Bi collision, at 3 GeV proton beam 
energy; results of the HSD model calculations 

the nucleus in the beam direction behaves very similar during first stage of 
the reaction as the average value of excitation energy. Starting from about 
10 fm/c, when proton reaches target nucleus, up to about 18 fm/c, the longi- 
tudinal momentum per nucleon of the residual nucleus has a maximal value, 
adequate to the value of momentum introduced into target nucleus by incom- 
ing proton. Then, the average value drops quickly in the range 20 - 30 fm/c 
of the duration time of the first stage of the reaction, because the momen- 
tum is carried out by nucleons escaping from the nucleus. At larger times 
the value varies only a little. 

Very similar behavior is observed for the average value of angular momentum 
of the residual nucleus during first stage of the reaction, presented in Fig. 
13.31 At about 10 fm/c, incoming proton, having in average a non zero impact 
parameter, introduces some angular momentum into the nucleus. At around 
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Figure 3.2: Time variation of the average value of longitudinal momentum 
per nucleon of the residual nucleus in p + Bi collision, at 3 GeV proton beam 
energy; results of the HSD model calculations 
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Figure 3.3: Time variation of the average value of angular momentum of the 
residual nucleus in p + Bi collision, at 3 GeV proton beam energy; results of 
the HSD model calculations 



18 fm/c and later the angular momentum lowers significantly due to nucle- 
ons escaping from the nucleus. This leakage of angular momentum stops at 
around 30 fm/c and the average value stabilizes. At later time, i.e. from 

about 40 fm/c, spurious slow increase of angular momentum is observed, 
which is unphysical and results from building-up inaccuracies of numerical 



Chapter 3. Specific models 



35 



calculations. That also indicates that the first stage calculations should be 
terminated at about 35 fm/c. 

Looking at the time evolution of the average mass number of the residual 
nucleus, presented in Fig. 13.41 quite different behavior is observed. It is seen, 
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Figure 3.4: Time variation of the average value of mass number of the residual 
nucleus in p + Bi collision, at 3 GeV proton beam energy; results of the HSD 
model calculations 

that at about 10 fm/c the projectile come into the target nucleus. As ex- 
pected, the mass number is increased by one. But then, starting from about 
18 fm/c, the mass number of the residual nucleus decreases monotonically 
with duration time of the first stage of reaction and does not stabilize at 
larger times, as it was in case of the other quantities. This indicates, that 
emission of particles in this model takes place all the time. Duration time of 
the cascade of the nucleon - nucleon collisions cannot be determined on the 
basis of behavior of the average mass number of the nucleus. Nevertheless, 
it accents the importance of assumption of a proper stopping time for the 
calculations, which has infiuence on a value of multiplicity of emitted parti- 
cles. It has a negligible meaning for energy distributions of ejectiles, since, 
as a time evolution of excitation energy shows, from about 18 fm/c, particles 
with lower and lower energies are emitted. 

The behavior of excitation energy, longitudinal momentum and angular mo- 
mentum as a function of time shows that choice of the stopping time is limited 
from two sides. The first stage calculations cannot be terminated too early, 
when interacting system is highly nonequilibrated. From the other side, also 
not too late, since for longer times numerical inaccuracies are increasing. 
Based on above dependences, it is concluded that the most reasonable dura- 
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tion time of transport models calculations of the first stage of proton induced 
reactions is equal to 35 fm/c. 

Behavior of the time dependences in case of reactions on other targets 
nuclei and other values of incident energy is displayed in Fig. 13.51 It is seen 
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Figure 3.5: Time variation of the average values of excitation energy per 
nucleon and angular momentum of residual nuclei in proton induced reactions 
on various targets, at different projectile energies; results of the HSD model 
calculations 



that stabilization of the average values in time depends both on incident en- 
ergy and mass of target. The heavier target is used, the later stabilization 
occurs. Similarly, the lower projectile energy, the longer time of first stage 
calculations is needed in order to reach equilibrium. Unfortunately, it cannot 
be assumed one maximal time for all systems, because of building up of the 
numerical inaccuracies (see the time dependences of angular momentum), 
what starts at different time for different systems. Nevertheless, the deduced 
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above stopping time equal to 35 fm/c will be treated in the frame of this work 
as an optimal time for each system. It is optimal for all used targets, since 
it has been established on the example of one of the heaviest targets. But 
taking into consideration dependence on incident energy, the following excep- 
tions must be included. Calculations of first stage of reactions at projectile 
energies greater than or equal to 3.0 GeV, on all targets will be terminated 
after 35 fm/c. At energies lower than 3.0 GeV, the HSD calculations will be 
stopped later, especially for heavier targets. 

The above considerations are based on the results of calculations averaged 
over all impact parameters. For particular reaction p+Au at Tp=2.5 GeV, 
for specific impact parameters, e.g. b = 1 fm, 4 fm, 6 fm, the following 
dependences: time evolution of average excitation energy (presented in Fig. 
13.61) . angular momentum (Fig. 13.71) and momentum in beam direction (Fig. 
13.81) are obtained. 
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Figure 3.6: Time variation of tlie average value of excitation energy per 
nucleon of the residual nucleus in p + Au collision at 2.5 GeV proton beam 
energy, the impact parameter b=6 fm; results of the HSD model calculations 

The example dependences calculated for specific impact parameters agree 
with results obtained with averaging over all impact parameters; choice of 
time of calculations of first stage of the reaction (^ 35 - 45 fm/c) is satisfac- 
tory. 




38 



3.5 Stopping time criteria 





80 










70 






p+ ^^''AuTp=2.5GeV 




60 




— X 








b = 4 fm 




50 








n 












40 








A 










—I 

V 


30 










20 










10 




















10 20 30 40 50 60 
time (fm/c) 

Figure 3.7: Time variation of the average value of angular momentum of the 
residual nucleus in p + Au collision at 2.5 GeV proton beam energy, the 
impact parameter b=4 fm; results of the HSD model calculations 




Figure 3.8: Time variation of the average value of longitudimal momentum 
of the residual nucleus in p + Au collision at 2.5 GeV proton beam energy, 
the impact parameter b—1 fm; results of the HSD model calculations 



Chapter 4 



Bulk properties of the first 
stage of proton induced 
reactions 

The following scenario of the first stage of proton - nucleus reactions takes 
place. High energy proton (Tp > 100 MeV) hits a target nucleus. The intra- 
nuclear cascade starts to develop: the incoming proton, on its way, collides 
with several target nucleons, transfers energy and momentum to them and 
may excite them into higher baryonic states. Depending on the position of 
the cascade nucleons inside the target, they either escape directly from the 
nucleus or collide secondarily with a few other nucleons transfering further 
the energy and momentum. Some of the nucleons may leave the nucleus. 
Development of the cascade can be seen by observing variations of spatial 
nucleon density with the reaction time. It indicates that proton induced 
reactions are low-invasive processes. 

In Figure 14.11 the HSD simulations of time evolution of nucleon density 
in central collisions of 2.5 GeV and 5.0 GeV proton with Au nucleus are 
presented. The density of nucleons with kinetic energies Ek > 50 MeV has 
been calculated, in order to see more clearly the development of the intra - 
nuclear cascade. Initially, before proton strikes the target nucleus, the max- 
imal kinetic energy of nucleons in the target is equal to the Fermi energy 
{Ep pa 35 MeV). As a consequence of interaction nucleons get more ener- 
getical. At the beginning of reaction, the incident proton is in distance of 
15 fm from the center of target nucleus. Projectile enters the target from 
the left, z is the axis in beam direction. One can see in the Fig. 14.11 that 
after 5 fm/c (Ifm/c ~ | ■ 10~^^ s), when proton has not entered the target 
yet, none of target nucleons has kinetic energy more than 50 MeV. Then, 
looking at the situation corresponding to density distribution sampled after 
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Figure 4.1: Time evolution of nucleon density of nucleons with kinetic ener- 
gies Ek > 50 MeV, in central collisions of 2.5 GeV (left column) and 5.0 GeV 
(right column) proton with Au nucleus; results of the HSD model calculations 
(projections on xz plane) 
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15 fm/c, one observes something like a wave going near the surface of target 
nucleus. Mainly nucleons from the most outer area of nucleus take part in 
the cascade. The center of target is not touched yet. Situation after 35 fm/c, 
it means after first fast stage of the reaction is shown on the bottom plots 
of the Fig. 14.11 It is seen that the center of the nucleon density distribution 
is now occupied by the cascade nucleons. But the density is not spatially 
uniform. The central area corresponds to the maximum of the distribution. 
Levels of constant nucleon density (in units nucleon/fm'^), with values in- 
creasing to the center of the circular shaped distribution are clearly visible. 
It means that in the residual nuclei, the density distribution of shape like of 
initial target nucleus is reproduced (the density of initial nucleus have the 
Wood-Saxon distribution form 03.131) ). Nevertheless, a small expansion of 
the nucleus is observed. The level corresponding to the spatial density of 
value less than about 0.002 nucleon/fm^ is associated with the free nucleons 
knocked out of the target. 

Distribution of the nucleon density of residual nuclei, in general, does not 
depend significantly on the centrality of proton - nucleus collision. This is 
illustrated in Fig. 14. 2[ where two-dimensional projections of the nucleon den- 
sity of residual nuclei, for different centralities of 5.0 GeV proton with Au 
nucleus collision are compared (i.e. for impact parameters of b=1.0, 2.5 and 
5.0 fm). However, on the level corresponding to the lowest spatial density, 
associated with emitted nucleons, some differences are observed. They indi- 
cate that more nucleons are emitted in central than in peripheral collisions. 

All the presented distributions of nucleon density show that in general, 
the incoming proton has caused only minor changes of density inside target 
nucleus. 

One of characteristics of created target nucleus is its maximal density. 
Modification of the density due to penetration of the nucleus by incoming 
proton is the characteristic feature of proton induced reactions. Such modi- 
fications are presented qualitatively in Fig. 14.31 as time evolution of a ratio 
of maximal nuclear density in case of p+Au reaction at 0.5 GeV, 1.0 GeV, 
2.5 GeV and 5.0 GeV of incident energy and the standard nuclear density 
Pmax- is evident from the Figure, that the incoming proton causes negli- 
gible modifications of nuclear density. Though, it is seen, that in case of low 
energy projectile the maximal density is first slightly increased and then it 
decreases. In case of higher energetical proton, first slight decrease and then 
an increase is observed. That could be explained by a fact that incoming 
low energy projectile almost stops inside target nucleus, causing increase of 
density, while higher energy projectile goes faster through the nucleus, push- 
ing nucleons away. Then, the situation is changed respectively, due to acting 
of mean field potential. Nevertheless, the deviations of the presented ratios 
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from the unity are of order of few percent, what proves that proton induced 
reactions are quite non-invasive processes. 

In order to check, how the kinetic energy is distributed inside the residual 
nuclei, the following test has been made. Based on p+Au reaction, at 2.5 
GeV proton beam energy, nucleons, which take part in the intra-nuclear cas- 
cade have been observed. It means, nucleons which gained a significant part 
of the energy of the projectile (i.e. of > 75 MeV, Ef, > 100 MeV). Most 
of the observed cases indicated homogeneous distribution of kinetic energy 
inside the residual nuclei. Partial heating up of the nuclei has been noticed 
in only about 1% of the cases, in peripheral collisions. Illustration of the 
exceptional situation is shown in Fig. 14. 4[ where, starting from the top, the 
nucleon density of residual nuclei of nucleons with kinetic energies Ei^ > 50 
MeV, Ek > 75 MeV and Ej, > 100 MeV, respectively are plotted. Looking at 
the plot in the bottom of Fig. 14.41 one can distinguish two differently excited 
parts. One part of nucleus is composed of nucleons with kinetic energies 
greater than 100 MeV, the second part - of nucleons with kinetic energies 
lower than 100 MeV. The highly excited group of nucleons, it is about 10 nu- 
cleons with an average total momentum equal to about 560(+/-80) MeV/c. 
The less excited part, it is till about 184 nucleons with an average total mo- 
mentum equal to only about 190(+/-60) MeV/c. 

Similar situation is observed, if looking at a bit lower energetical nucleons, 
the central plot of Fig. 14.41 nucleons with kinetic energies greater or lower 
than 75 MeV. In this case, the highly excited part consists of about 20 nucle- 
ons with an average total momentum equal to about 470(+/-90) MeV/c, and 
less excited part composed of 174 nucleons with an average total momentum 
equal to about 190(+/-60) MeV/c. 

One can conclude, that the HSD simulations predict something like forma- 
tion of two excited sources of evidently unequal masses. The smaller source 
consists of relatively few nucleons (up to ~ 20) and is rather fast, (3 ~ 0.025 c. 
The larger source is built of ~ 170 - 180 nucleons and has velocity (3 ~ 0.0012 
c. Similar observation has been drawn from phenomenological analysis of ex- 
perimental data presented in 
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Figure 4.2: Nucleon density of residual nuclei of nucleons with kinetic en- 
ergies Ek > 50 MeV, in central (starting from the top, b=lfm), midcentral 
(b=2.5fni) and peripheral (b=5fni) collisions of 5.0 GeV proton with Au nu- 
cleus; results of the HSD model calculations (the density projections on xz 
(left column ) and yz (right column ) planes ) 
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Figure 4.3: Modifications of maximal nucleoli density in the case of p+Au 
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results of the HSD model calculations 
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Figure 4.4: Nucleon density of residual nuclei of nucleons with kinetic energies 
Ek > 50 MeV (starting from the top), > 75 MeV and Ek > 100 MeV, 
respectively, in peripheral (b—5fm) collisions of 2.5 GeV proton with Au 
nucleus; results of the HSD model calculations 



Chapter 5 



Properties of residual nuclei 
after the first stage of proton - 
nucleus reactions 

As it is mentioned above in this work, as result of the first stage of proton 
- nucleus reaction, beside emitted particles, an excited nucleus remains. It 
differs from the initial target, in average, by only a few nucleons in mass 
number. Properties of the residual nucleus (i.e. mass (Ar), charge {Zr), ex- 
citation energy three - momentum (pjj) and angular momentum (Lpi)) 
are evaluated, in the frame of the HSD model, by exploring conservation 
laws, according to formulas 13.231 Results of calculations for reactions on 
various target nuclei, at different values of incident energy in range from 0.1 
GeV to about 10 GeV, are discussed below. 

Let's look first at one - dimensional distributions of the quantities evolving 
with projectile energy and mass of target. Histogrammed properties for ex- 
emplary reactions of 1.9 GeV proton on light (^^Al), heavy (^^^Au) and two 
intermediate mass (^^Ni and ^'''^Ag) targets are displayed in Fig. 15. 1[ Distri- 
butions for proton induced reaction on example Al target, at several values of 
proton beam energy are shown in Fig. 15. 2[ The average values and standard 
deviations for the presented distributions are collected in Tables 15.11 and 15.21 
The Figures 15.11 and 15.21 and values collected in the Tables 15.11 and 15.21 
indicate, that all of the distributions differ significantly with mass of target, 
but behave similarly for varied projectile energies. 
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Table 5.1: The average values and standard deviations for the one - di- 
mensional distributions of properties of residual nuclei from 1.9 GeV proton 
induced reactions on several targets 



reaction 


P+Al 


P+Ni 


P+Ag 


p+Au 


<A>B. 


23.96±2.25 


54.71±2.65 


102.24±3.55 


192.39±3.69 


< Target Mass loss > 


3.036±2.25 


3.29±2.65 


4.76±3.55 


4.61±3.69 


<Z>R 


11.64±1.42 


26.53±1.64 


45.05±1.98 


77.34±1.92 


< Target Charge loss > 


1.36±1.42 


1.47±1.64 


1.96±1.98 


1.66±1.92 


< E* >R [MeV] 


74.56±81.77 


122.84±120.41 


144.16±125.65 


222.16±171.24 


< E*/N >R [MeV/N] 


3.27±3.63 


2.30±2.29 


1.43±1.27 


1.16±0.91 


<L>R [h] 


4.77±3.25 


7.41±3.59 


11.76±5.50 


17.14±8.08 


<p, >R [GeV/c] 


0.18±0.24 


0.22±0.27 


0.29±0.35 


0.34±0.38 


< >R [GeV/c] 


-0.0011±0.28 


-0.0014±0.31 


-0.0058±0.36 


-0.0045±0.37 



Table 5.2: The average values and standard deviations for the one - di- 
mensional distributions of properties of residual nuclei from proton induced 
reactions on Al taxget, at several incident energies 



impact energy 


1.0 GeV 


2.0 GeV 


3.0 GeV 


4.0 GeV 


<A>R 


23.51±2.17 


24.032±2.20 


24.23±1.99 


24.12±2.085 


<Z>R 


11.28±1.35 


11.66±1.37 


11.78±1.31 


11.71±1.33 


< E* >R [MeV] 


65.76±54.067 


73.46±84.27 


70.31±88.081 


74.50±98.023 


< E*/N >R [MeV/N] 


2.94±2.61 


3.22±3.75 


3.026±3.84 


3.25±4.35 


<L>R [h] 


4.75±2.65 


4.71±3.072 


4.85±4.48 


4.99±5.23 


< p, >R [GeV/c] 


0.18±0.24 


0.17±0.23 


0.16±0.24 


0.16±0.25 


< >R [GeV/c] 


-0.0021±0.27 


-0.00088±0.27 


-0.0030±0.27 


0.0033±0.27 


impact energy 


5.0 GeV 


6.0 GeV 


7.0 GeV 


8.0 GeV 


<A>R 


24.027±2.19 


23.94±2.26 


23.86±2.35 


23.75±2.41 


<Z>R 


11.66±1.39 


11.603±1.41 


11.58±1.46 


11.54±1.47 


< E* >R [MeV] 


74.61±90.60 


76.48±98.98 


79.64±129.89 


81.92±137.73 


< E*/N >R [MeV/N] 


3.29±4.13 


3.41±4.55 


3.58±5.8 


3.71±6.29 


<L>R [h] 


5.034±5.66 


5.48±8.25 


5.15±7.32 


6.035±11.044 


< p, >R [GeV/c] 


0.16±0.24 


0.16±0.25 


0.17±0.27 


0.18±0.28 


< Px >R [GeV/c] 


-0.000024±0.27 


0.0053±0.28 


0.0026±0.28 


-0.0051±0.29 
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Figure 5.1: One - dimensional distributions of properties of residual nuclei 
after first stage of 1.9 GeV proton induced spallation reactions on '^'^Al, ^^Ni, 
^°''Ag and ^^''Au targets; results of the HSD model calculations 
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Figure 5.2: One - dimensional distributions of properties of residual nuclei 
after first stage of proton induced spallation reactions on '^'^ Al target, at few 
values of projectile energy; results of the HSD model calculations 
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Reason for these behaviors is connected with finite size of nuclei. Average 
values of excitation energy (E^) and angular momentum (Lp/) of residual 
nucleus and width of the distributions increase with mass of target nuclei, 
at fixed projectile energy. For a chosen target, the distributions of very 
slightly depend on incident energy. Distributions of almost do not depend 
on incident energy. This can be explained in an intuitive way by analyzing 
how the energy and momentum are deposited into a target nucleus. The 
E"^ left in the residual nucleus after the first stage of reaction is in average 
only a small fraction of the energy, which has been introduced into a target 
nucleus by incoming proton. Most of the energy is carried away by emitted 
particles (example energy balance is presented in Sec. 18. ip . The number of 
ejected particles should be treated as a signature how many nucleons were 
involved in the intra - nuclear cascade. The more nucleon - nucleon collisions 
inside target nuclei, the larger part of the energy and angular momentum 
introduced by a projectile is deposited into the residual nucleus. E.g. due to 
a collision of the incident proton with a nucleon inside a target, a resonance, 
mainly Delta resonance could be excited. The Delta resonance has a short 
time of life, i.e. about 1-2 fm/c. It decays into nucleon and meson (mainly 
pion) or collides with other nucleon before the decay. As result, energy is 
accumulated inside target nucleus. The heavier target, the more collisions can 
occur, because the cascade of collisions can last for a longer time. For light 
targets, the excitation energy (E^) is respectively small, because a number 
of collisions made by a projectile is low. In case of the excitation energy per 
nucleon (E^/N) the situation is opposite, i.e. E^/N is a decreasing function 
of mass of target. This is just because in heavier targets, the amount of 
energy is distributed among more nucleons. 

Similarly, distribution of angular momentum of residual nuclei almost 
does not depend on value of incident energy, but strongly depends on mass 
of target. The heavier target, the more angular momentum is deposited 
inside. This is because in the HSD model, quasi - classical approximation 
for calculations of angular momentum is used; at fixed impact energy, the 
maximal angular momentum of projectile is proportional to radius of target 
nuclei. Therefore, in case of heavier targets, more angular momentum is 
introduced by a projectile. 

Distributions of momentum in beam direction (z-component) and the 
momentum perpendicular to the beam direction (x-component) indicate that 
residual nuclei are moving according to initial direction of projectile. It is 
evidently seen (both from the Fig. 15. H and Table [5?T1) . that the average values 
of momentum in x-direction of residual nuclei are consistent with zero. 
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Table 5.3: The average values and standard deviations for distributions of 
velocity of residual nuclei (< >r [c]) from proton induced reactions on 
several targets 



<V^>R [c] 


p+Al 


p+Ni 


p+Ag 


p+Au 


Tp=1.2 GeV 


0.0082±0.011 


0.0046±0.0056 


0.0036±0.0036 


0.0035±0.0027 


Tp=1.9 GeV 


0.0078±0.011 


0.0044±0.0054 


0.0029±0.0033 


0.0025±0.0024 


Tp=2.5 GeV 


0.0069±0.010 


0.0036±0.0052 


0.0024±0.0032 


0.0019±0.0021 



Similar information about the reaction could be extracted from distribu- 
tions of velocity in beam direction of residual nuclei. The distributions differ 
significantly with mass of target, but behave similarly for varied projectile 
energies, as shown in Fig. 15.31 The average values and standard deviations 
of the values are presented in Table 15. 3[ It is seen that the lighter target, 
the broader is the distribution. These are asymmetric, peaked in forward 
direction, with tails in backward direction. This indicates, that in average, 
the residual nuclei tend to move forward, according to the initial direction of 
projectile. 

In order to have a better view of the dependences, average values of the 
quantities and standard deviation of the average values, have been expressed 
as a function of incident energy and mass of target and presented in Figures 
OlandOl 

Average values and standard deviations of mass number and charge of 
residual nuclei remaining after first stage of p+Al reaction, in function of 
projectile energy are presented in Figures 15.61 and 15.71 It is interesting, that 
both the mass number and the charge, in average, increase up to 3.0 GeV 
of incident energy. Then decreasing is observed. The ratio of average values 
of the mass number and charge is constant as function of projectile energy, 
as shown in Fig. 15. 8[ Additionally, it exceeded the ratio for initial target 
nucleus. 

Average values and standard deviations of momentum in beam direc- 
tion of residual nuclei after first stage of proton induced reaction on several 
targets, at various incident energies are presented in Fig. 15.91 One can con- 
clude, that the heavier the target, the more momentum in beam direction is 
deposited inside. What more, the maximum in mass dependence on projec- 
tile kinetic energy (at the incident energy equal to 3.0 GeV) for given target, 
see Fig. 15.61 corresponds to the minimum of dependence of momentum in 
beam direction of residual nuclei in function of projectile energy. 
Behavior of standard deviations as function of projectile kinetic energy is 
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noteworthy. It is observed, that the higher incident energy, the broader 
are distributions of excitation energy, momentum, angular momentum, mass 
number, charge of residual nuclei. This can be qualitatively explained taking 
into account results of simulations presented in the Fig. 15.101 (i.e. depen- 
dence of multiplicity and energy of emitted nucleons during the first stage 
of the reactions on incident energy). The behavior of multiplicity of emitted 
nucleons as function of projectile energy perfectly agrees with the behav- 
ior of mass and charge of residual nuclei. It means, first, up to about 3.0 
GeV of incident energy, number of emitted nucleons decreases, so that mass 
of residuum increases. With further increase of impact energy, more and 
more nucleons are emitted, so that mass of the residuum decreases. Obvi- 
ously, average kinetic energy carried by emitted nucleons is a monotonically 
increasing function of projectile energy in the whole considered range of pro- 
jectile energy, as shown in Fig. 15.101 Protons carried out in average more 
kinetic energy than neutrons, since more protons than neutrons are emitted 
during first stage of reaction. The p/n ratio is larger than unity, because of 
the extra proton correspondig to the incident particle (see Fig. 15.101) . One 
can conclude that information of type of projectile is somehow remembered 
by the system during the first stage of reaction (this topic is undertaken also 
in Sec. E2). 

Looking at behaviors of average values of properties of residual nuclei, 
especially their momentum in beam direction, mass number and multiplicity 
of emitted nucleons during first stage of reaction as function of incident en- 
ergy, it is seen that the average values first continuously decrease (or increase, 
respectively) with the incident energy, up to about 3 GeV and then consis- 
tently, the opposite tendency is observed. The precise explanation of such 
behavior is difficult, it must be a result of a few simultaneous effects. New- 
ertheless, considering elementary processes (implemented to the model, based 
on experimental data), one can suggest the following qualitative explanation. 
Such behaviors should be connected with contribution of different resonances 
excitations, which depends on the energy introduced by a projectile into the 
target nucleus. At low projectile energy (about 1.0 GeV) nucleons are emit- 
ted mainly due to elastic nucleon - nucleon scattering. More probable than 
at higher energies is also, that incoming proton stops inside target nuclei. At 
higher incident energies possibility of excitation of resonances (mainly, the 
most important here A resonance) increases. A bit less, but more energetic 
nucleons are emitted. In result, the average multiplicity of emitted nucleons 
decreases, so mass of residual nuclei increases and the average longitudinal 
momentum decreases. At higher incident energy (equal to about 3.0 GeV), 
from one side the probability of A resonance excitation becomes smaller, and 
from the other side probability of excitation of higher resonance is still too 
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low. This corresponds to the minimum (or maximum) of the incident en- 
ergy dependences of the considering quantities. Then, with further increase 
of projectile energy, probability of excitation of other resonances increases. 
This causes, that more energy and momentum is deposited into residual nu- 
clei. So, the average longitudinal momentum of residual nuclei increases. 
Simultaneously, more particles could be produced, since the ratio of inelastic 
to elastic collision increases. The multiplicity of emitted nucleons increases, 
so the average mass of remnants decreases. 

Behavior of multiplicity and average kinetic energy of emitted nucleons, as 
function of target mass number, at an example impact energy is shown in 
Fig. 15.111 One can see, that the heavier target, the smaller differences in 
number of emitted protons and neutrons. For the heaviest presented target, 
the p/n ratio is equal to the unity. One can conclude that in heavy targets, 
the information about a type of projectile is lost due to possible large num- 
ber of mutual reactions. Total number of emitted nucleons increases with 
mass of target nucleus and stabilizes for heavier targets (with mass numbers 
greater than 60). In heavier targets, the cascade of intra- nuclear collisions 
can last longer, so that emitted nucleons carries in average less kinetic en- 
ergy. The nucleon, on its way, going outside the nucleus, has to penetrate 
through nuclear skin of the nucleus, with high probability undergoes further 
collisions. As result, more particles with suppressed energies are emitted. 

More informations about the excited nuclei is provided by two - dimen- 
sional correlation plots of their properties. Example dependences of excita- 
tion energy, mass of residual nuclei and their velocity in the beam direction 
are shown in Fig. 15.121 (velocity versus excitation energy per nucleon and 
velocity versus residual mass). It is meaningful, that velocity in beam direc- 
tion of residual nuclei increases with increasing excitation energy deposited 
inside the residual nuclei. This indicates, that the energy is comprised also in 
onward movement of residual nuclei. Simultaneously, the velocity increases 
with decreasing mass of residual nuclei, what is equivalent to increase of 
number of ejected particles. One can conclude that also excitation energy of 
residual nuclei should increase with decrease of mass of the nuclei. In Figures 
15.131 and 15.141 regression function of excitation energy on target mass loss 
and regression function of angular momentum on target mass loss are pre- 
sented (for residual nuclei formed in p-l-Al reaction, at a few impact energies 
and proton induced reaction on several targets at example 4.0 GeV incident 
energy, respectively). The dependences are roughly linear functions of the 
target mass loss. 

It is seen, that the heavier remnants (for the same initial target), the 
lower average excitation energy is deposited inside. The average excitation 
energy increases with decreasing of the remnant mass. Based on the Fig. 
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15.14^ one can conclude that particularly the dependence on mass of target is 
quite pronounced. It is because, as it is already mentioned above, in heavier 
targets, where the cascade lasts longer, more energy can be accumulated. In 
the Figures [5.131 and 15. 141 it is also seen, that in the cases of the lowest mass 
of residual nuclei, (i.e. the target mass loss higher than about 7 nucleons), 
the increase of excitation energy and angular momentum is not only less 
pronounced, but even a slight decreases is observed. In such cases, i.e. more 
ejected particles, the intra - nuclear cascade is lasting for longer time, so that 
particles with rather moderated energies are emitted. 

Example two - dimensional distributions of mass and excitation energy 
of residual nuclei formed in p+Al reaction, at several projectile energies are 
displayed in Fig. 15.151 Again, one observes that the higher excitation energy 
of residual nuclei, the smaller their mass. 

Dependence of excitation energy on angular momentum of residual nuclei 
is presented in Fig. 15.161 The five lines (in upper part of the Fig. 15.161) 
correspond to the Yrast lines (definition is given in Chapter [Tj) calculated 
according to the formula: E* = ^Lr{Lr + 1) with different values of the 
moment of inertia (J), where: I is the moment of inertia of rigid body (/ = 
|Mi?^, where M and R stand for mass and radius of nucleus, respectively), 
//50 corresponds to the hydrodynamical limit of nucleus (meaning that only 
a part of nucleus behaves like the rigid body, the other part - in uncorrelated 
way, see Ref. [52]). One observes that angular momentum increases with 
increase of the average value of excitaton energy. It is evident, that the 
dependence is not inconsistent with the Yrast line. 

Physicists, who are studing fragmentation have found that this process 
occurs at some particular value of projectile energy (Tp) (based on (p,^Be) 
reaction, for various targets at different Tp, see [S3]). Fragmentation is a 
phase transition, what means, that a proper amount of energy needs to be 
pumped into a system. It follows, that fragmentation occurs at some partic- 
ular value of {E* /N)fi, which experimentally has been found to be at least 
equal to about 5 MeV/N [Ml [55]. Using the condition of 5 MeV/N to the 
HSD model results, presented in Fig. 15. 17^ it is seen that fragmentation is 
a very low probable process in proton induced reactions. For example, in 
p+Al reaction at Tp=5.0 GeV, fragmentation of excited nucleus may occur 
in about 17 % of cases (i.e. the part of spectrum for {E* /N)fi. > 5 MeV/N 
is 17 % of the whole {E* /N)ji spectrum), whereas in p+Au reaction at the 
same incident energy, only in about 5 % of cases. These are cases from the 
tail of {E*/N)ji distribution. Based in the Fig. 15.171 one can conclude that 
fragmentation process is more probable in reactions on light than on heavy 
targets, and at higher incident energies. 
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One can suppose that fragmentation of only a part of nucleus could be 
more probable. This requires inhomogeneous excitation of residual nucleus. 
It means, not all of nucleons inside the nucleus are involved in carrying of 
the excitation, i.e. for different target, {E* /N)ii should be calculated not 
per each nucleon, but only per some amount of nucleons, which take active 
part in the cascade. However based on observations of evolution of nucleon 
density and distribution of kinetic energy inside excited nuclei, it is seen that 
such partial heating up of nucleus is very low probable. It occurs only in 
about 1 % of discussed in Chapter HI 

Additionally, experimentalists, who are interested in fragmentation, look 
for it in central collisions, accompanied with high multiplicity of emitted par- 
ticles, as it is shown qualitatively by results of the HSD model calculations, 
presented in Fig. 14. 2[ Nevertheless, in the HSD simulations, the very low 
probable partial heating up of nucleus has been found in peripheral collisions 
(see Fig. I44|l . 

Investigated here reactions are very low invasive processes, involving only 
a few nucleons, causing only minor fluctuations of nucleon density. So, it is 
evident that there is a negligible probability for fragmentation in proton - 
nucleus interacting system, in energy range considered in the frame of this 
work. Spallation is evidently dominant process. 
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Figure 5.3: One - dimensional distributions of velocity in beam direction of 
residual nuclei at the end of first stage of proton induced spallation reactions 
on various targets, at few example values of projectile energy; results of the 
HSD model calculations 
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Figure 5.7: Standard deviation of average values of mass number (left) and 
charge (right) of residual nuclei from proton induced reactions on Al target, 
in function of incident energy; results of the HSD model calculations 
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Figure 5.8: Ratio of average values of mass number and charge of residual 
nuclei from proton induced reactions on Al target, in function of incident 
energy; results of the HSD model calculations 
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Figure 5.9: Average values (left) and standard deviation (right) of the average 
values of momentum in beam direction of residual nuclei from proton induced 
reactions on several targets, in function of incident energy; results of the HSD 
model calculations 
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Figure 5.10: Multiplicity (left) and average values of kinetic energy per indi- 
vidual nucleon (right ) of nucleons emitted during first stage of proton induced 
reaction on Al target, as function of incident energy; results of the HSD model 
calculations (error bars indicate values of standard deviation of the average 
values, divided by a square root of number of events) 
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Figure 5.11: Multiplicity (left) and average values of kinetic energy per in- 
dividual nucleon (right) of nucleons emitted during first stage of 4.0 GeV 
proton induced reactions on a few targets, as function of target mass num- 
ber; results of the HSD model calculations (error bars indicate values of 
standard deviation of the average values, divided by a square root of number 
of events) 
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Figure 5.12: Two - dimensional distributions of velocity in beam direction 
of residual nuclei and (a) excitation energy per nucleon of residual nucleus 
and (b) mass of residual nuclei, for exemplary reaction p+Al at 1.2 GeV 
of incident energy; results of the HSD model calculations (the most central 
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Figure 5.13: Regression function of excitation energy (left) and angular mo- 
mentum (right) of residual nuclei on target mass loss, for example p+Al 
reaction, at several values of incident energy; results of the HSD model cal- 
culations (error bars indicate values of standard deviation of average values of 
excitation energy and angular momentum, respectively, for adequate target 
mass loss, divided by a square root of number of events) 
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Figure 5.14: Regression function of excitation energy (left) and angular mo- 
mentum (right) of residual nuclei on target mass loss, for proton induced 
reactions on Al, Ca, Ni, Ag and Au nuclei, at an example value 4.0 GeV of 
incident energy; results of the HSD model calculations (error bars indicate 
values of standard deviation of average values of excitation energy and an- 
gular momentum for adequate target mass loss, divided by a square root of 
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Figure 5.15: Contour plots of the two-dimensional distributions of excitation 
energy and mass of residual nuclei at the end of first stage of proton induced 
reaction on Al target, at several values of incident energy; results of the HSD 
model calculations (the most central contour line corresponds to a maximum 
of the distribution, the more and more outer contours correspond to the 
decrease of a yield) 
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Figure 5.17: Distributions of excitation energy per nucleon of residual nuclei 
from an example proton induced reaction on Al and Au target, at 2.5 GeV 
and 5.0 GeV of incident energy; results of the HSD model calculations 
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5.1 Parametrization 

It is impossible to measure or even to calculate the reaction of proton of all 
energies on all target nuclei. One comes to idea to find a global parametriza- 
tion for distributions of excitation energy versus mass of residual nuclei for 
some chosen sample of targets and energies, involving finally dependency on 
the target mass and projectile energy, in order to have the possibility to in- 
terpolate for other energies and targets. 

Two-dimensional distributions of excitation energy and mass of residual nu- 
clei remaining after first stage of proton induced reaction, on a few targets, 
at several projectile energies have been parametrized by a sum of at least 
three two - dimensional Gaussians: 

f{AR,E*)= ai■exp{-b,{AR-Ano^Y-Ci{E*-E*f) (5.1) 

i=0,l,2 

where: aj stands for heights of the Gaussians, bi and Cj - define widths in 
and E* directions, likewise Ajioi and Eq^ - the mean position of the 
Gaussians in Ar and E* directions, respectively. Example results of such 
parametrization are presented in Figures 15.181 and 15.191 It is seen, that the 
distributions, calculated with the HSD model, are very well described by a 
sum of Gaussians fl5.1l) . 

Values of parameters of the Gaussians evolve rather smoothly with inci- 
dent energy and mass of target nuclei. Example results of parametrization 
of such dependences are presented in Figures 15.201 and 15.211 

Using the parametrizations, two-dimensional distributions of excitation 
energy versus mass of residual nuclei could be obtained for any value of pro- 
jectile energy and target nuclei from the chosen ranges, for the considered 
cases. In order to verify the parametrization, interpolated distributions have 
been compared to results of the HSD calculations. Comparisons of distri- 
butions for example p+Al reaction at 2.5 GeV beam energy (interpolation 
of parameters in function of projectile energy, see Fig. I5.20p and p+Zr re- 
action at 2.0 GeV beam energy (interpolation of parameters in function of 
mass of target nuclei, see Fig. I5.2ip are presented in Fig. I5.22[ Obtained 
results are satisfactory. This indicates possibility of realization of such a 
parametrization for other target nuclei and other incident energies. 
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Figure 5.18: Contour plots of the two-dimensional distributions of excitation 
energy and mass of residual nuclei in p+Al reactions, at several values of 
beam energy; jagged lines correspond to results of the HSD model calcu- 
lations, smooth lines define parametrization (the most central contour line 
corresponds to a maximum of the distribution, the more and more outer 
contours correspond to the decrease of a yield) 
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Figure 5.19: Contour plots of the two-dimensional distributions of excita- 
tion energy and mass of residual nuclei in 2.0 GeV proton induced reaction 
on different targets; jagged lines correspond to results of the HSD model 
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Figure 5.20: Parametrization of values of the Gaussians parameters as func- 
tion of incident energy for the p+Al reaction; e.g. Aptoi = xl, Eq^ = yl, see 
formula 15. Jl (x^ means chisquare per degree of freedom ) 
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Figure 5.21: Parametrization of values of the Gaussians parameters as func- 
tion of target mass number for the 2.0 GeV proton induced reaction; e.g. 
Arqi = xl, = yl, see formula l5l] (x^ means chisquare per degree of 
freedom ) 
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Figure 5.22: Comparison of interpolated distributions and results of the HSD 
calculations (smooth and jagged lines) for example p+Al reaction at 2.5 GeV 
beam energy and p+Zr reaction at 2.0 GeV beam energy 
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Pions are the most abundantly produced mesons during proton induced re- 
actions on atomic nuclei. Creation of pions is very important process in 
the reactions due to the fact that they carry away significant part of four- 
momentum introduced by incoming particle during the first stage of the re- 
action. This is because pion mass (~ 140 MeV) is small enough to allow for 
copious production of pions in nucleon-nucleon collisions even at low beam 
energies, and large enough to contribute significantly to the energy and mo- 
mentum transfer. 

It should be emphasized that in spite of the fact that momentum carried 
away by pions is several times smaller than appropriate momentum of nu- 
cleons (see e.g. the energy balance presented in Sec. 18. ip . the contribution 
of pions to the momentum and energy transfer may be comparable to that 
of nucleons. This is because pions, due to the low production threshold can 
be produced and reabsorbed several times during the reaction (threshold for 
charged pion production in free nucleon - nucleon collision is equal to about 
289 MeV; inside nucleus, it can be decreased even up to about 141 MeV 
due to the Fermi motion). Thus they influence the whole dynamical evolution 
of the reaction. Proper treatment of pion production and absorption may be 
crucial for realistic estimation of excitation energy and other properties of 
residual nuclei after the fast cascade of nucleon - nucleon collisions. 
Calculations of the first stage of the reaction are sufficient in order to receive 
realistic pion spectra, since all pions are produced only in violent nucleon - 
nucleon collisions, where the locally available amount of four-momentum is 
large enough. 

One can distinguish two components in pion angular and energy distribu- 
tions: high energy anisotropic part, i.e. pions emitted dominantly in forward 
direction and isotropic part consisted of low energy thermal pions emitted in 
the whole angular range (0° - 180°). 
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Pions emitted from the first chance NN colhsions have high energy, whereas 
those from secondary colhsions are less energetic. Building up of pion spectra 
in time is presented in Figures 16.11 and 16.21 Typical momentum distributions 
of pions emitted from 1.0 GeV proton induced reaction on Al target are 
shown; Pz corresponds to the momentum in beam direction. The time evo- 
lution of pion distributions is depicted more precisely in the Fig. 16.21 where 
double differential kinetic energy spectra of pions produced during the ex- 
ample reactions on light (Al) and heavy (Au) target, at low (1.0 GeV) and 
higher (2.5 GeV) projectile energy are displayed. It is seen, that in case 
of 1.0 GeV proton induced reactions, high energy anisotropic part of pion 
distribution is formed first. Then, the isotropic part is built up. Formation 
of the low energy part of pion spectra takes longer time, since it is result of 
several acts of absorption and emission of pions. In the Fig. 16.21 all pro- 
duced pions at the chosen reaction time are plotted, at further time some of 
them can be absorbed or absorbed and re-emitted in other direction, so the 
number of presented pions can be lower than at ealier time. After 35 fm/c, 
all presented pions are emitted. The number of produced pions is higher in 
case of 2.5 GeV proton induced reactions. Pion distributions are evolving 
with reaction time only slightly. At fixed energy of incident protons, pions 
produced on heavy target are more abundant than those produced on light 
target and low energy pions are dominant. Pion on its way can be several 
times absorbed and emitted (its mean free path is equal to 1.0 - 1.5 fm), 
its initial direction could be changed and its energy is suppressed. So, it is 
much more difficult for pion to be emitted in forward direction in case of 
heavy target, where it has to pass a longer distance inside nuclear matter. 
Therefore, high energy pions emitted in forward direction dominate in case 
of light target. Quantitative comparison of yield of the produced pions is 
displayed in Table 16. 1[ It is noteworthy, that almost all of the pions have 
been produced due to Delta resonance decay; only a very few pions have been 
produced due to N(1440) decay. 

Spectra of pions change significantly with projectile energy. It is seen in 
Fig. 16.31 where examplary distributions of transversal (Pt) versus longitudi- 
nal {Pz) momentum of neutral pions emitted during proton induced reaction 
on Ni target are plotted for various incident energies. These are results of 
HSD model calculations, where Pz corresponds to the pion momentum in 



It can be seen, how the evolution of pion spectra looks like with increase 
of incident energy. At low projectile energies pions have not enough energy 
to pass through a target nucleus, so they cannot be emitted in forward di- 
rection. In this case, only a small part of target is accessible for pions, the 



beam direction and Pt is defined as 
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Table 6.1: Time evolution of multiplicity of pions from p+Al and p+Au 
reaction at 1.0 GeV and 2.5 GeV proton beam energy, results of the HSD 
model calculations 





p(1.0 GeV)+Al 


p(2.5 GeV)+Al 




TT 




7r+ 


TT 


ttO 


7r+ 


15 fm/c: 


0.067 


0.18 


0.29 


0.38 


0.51 


0.50 


20 fm/c: 


0.071 


0.18 


0.29 


0.37 


0.47 


0.47 


35 fm/c: 


0.069 


0.17 


0.26 


0.37 


0.47 


0.46 




p(1.0 GeV)+Au 


p(2.5 GeV)+Au 




TT 




7r+ 


TT 


tt" 


7r+ 


15 fm/c: 


0.010 


0.045 


0.060 


0.35 


0.46 


0.41 


20 fm/c: 


0.085 


0.19 


0.26 


0.53 


0.60 


0.52 


35 fm/c: 


0.11 


0.17 


0.23 


0.49 


0.53 


0.46 



rest of nucleus corresponds to the so-called nuclear shadow. Higher pion 
momenta are associated with increasing incident energy. With increase of 
the incident energy, the anisotropic part of distribution increases. It means, 
more and more pions have chance to be emitted in forward directions. The 
isotropic, low - momentum part of pion spectra at high projectile energy cor- 
responds to thermal pions, which result from sequence of several absorptions 
and emissions. One should also notice that at high projectile energies (in the 
Fig. 16.31 - higher or equal to 2.0 GeV), the isotropic part of pion momentum 
spectra is symmetric around Pz = 0. At low projectile energies (in the Fig. 
16.31 - less than 2.0 GeV), the isotropic part of spectra, i.e. the most central 
contour line corresponding to a maximum of distribution is not symmetric 
around Pz = 0, but shifted to negative values of Pz', to large extent it is 
due to mentioned nuclear shadow. The isotropic part of pion spectra that 
consists of thermal pions has the Maxwell distribution. This is illustrated by 
Fig. 16. 4[ where the experimental distribution of negative pions, measured by 
Cochran et al. [57J at backward direction (150°), is presented. The following 
form of Maxwell function has been used for the parametrization: 
f{E^-) = Ao ■ {{E^~ + m^-f - ml.) ■ e"^-/^" 

where: E^^- and m-Tr- are the kinetic energy of the pions and their mass, 
respectively. Aq and Tq are fitted parameters. The following values of 
the parameters are obtained: Aq = 11476.1 ± 364.5 [GeV~^] and Tq = 
0.024 ± 0.00039 [GeV]. 

It is interesting to see, how the average yield of pions produced in pro- 
ton induced spallation reactions changes with projectile energy and mass of 
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Figure 6.1: Momentum distributions of pions emitted during p+Al reaction, 
at 1.0 GeV beam energy; Pz corresponds to momentum in beam direction; 
results of tlie HSD model calculations (scatter plot) 
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Figure 6.2: Double differential kinetic energy spectra of all kind pions pro- 
duced during p+Al and p+Au reactions, at 1.0 GeV and 2.5 GeV proton 
beam energy; results of tlie HSD model calculations 



target. The HSD model calculations predict that multiplicity of pions, in 
average, increases both with incident energy and mass of target. One can 
expect that due to charge conservation, the relative yield of produced pions 
should fulfill the relation: tt^ > > . The incident energy dependence 
(presented in Fig. 16.51) shows, that positive and neutral pion production is 
favoured, as expected. For heavier targets, where there are more collisions. 



Chapter 6. Pion spectra 



77 



multiplicity of positive pions decreases and for the heaviest target it is even 
lower than multiplicity of neutral and negative pions, as displayed in Fig. 
I6.6[ The target mass dependence shows, that reabsorptions and secondary 
collisions play an important role, especially for heavy targets. 
Looking at kinetic energy carried by individual emitted pions, it is seen (Fig. 
16.51 right) that all kinds of pions produced in proton induced reactions on 
a chosen target carry, in average, the same amount of kinetic energy. The 
higher energy of the incident proton, the bigger amount of kinetic energy car- 
ried out by pions. In reactions on heavier targets (Fig. 16. 6p . emitted pions 
carry away, in average, less kinetic energy. In case of the heaviest targets, 
the negative pions take the least kinetic energy of all pions. The suppression 
of kinetic energy of emitted pions from heavier target is an expected result of 
numerous possible reabsorptions and emissions of pions, which become less 
and less energetic. 

Inclusive differential kinetic energy spectra of pions emitted at different 
angles, calculated with the HSD model, have been confronted with available 
experimental data , as displayed in Fig. 16. 7[ for an example distributions 
of negative pions from p+Pb reaction, at 0.73 GeV of incident energy. One 
finds, that the calculations overestimate the experimental data. It is evi- 
dent, that the pion dynamics in the HSD model needs to be improved, i.e. 
the number of produced pions should be reduced. Pions are created mainly 
due to decay of Delta(1232) resonances. It means, pions production and ab- 
sorption proceeds mainly through the NN ^ NA and A ^ tcN reactions 
(see description of the model in Chapter [3]). Therefore, modifications of the 
cross section for the absorption of Delta resonance in reaction A^A NN 
and the Delta lifetime should influence the amount of produced pions. In 
order to estimate the influence, the following corrections have been tested. 
First, the Delta resonance lifetime has been enlarged. Comparison of the pion 
spectra calculated with HSD model, with enlarged Delta resonance lifetime 
and the experimental data [57] is presented in Fig. 16.81 (the calculations are 
performed with 30 times longer lifetime of Delta resonance). It is seen, that 
this correction reduces number of pions emitted mainly in backward direc- 
tions. This is due to the fact, that longer lifetime of Delta resonance implies 
smaller number of acts of absorption and emission of pions. If the cross sec- 
tion for the absorption of Delta resonance in reaction A^A — >■ NN is enlarged 
(by factor 20), the number of pions emitted in forward directions is reduced, 
as it is shown in Fig. 16.91 Implementation of both of the corrections reduces 
the amount of produced pions both in forward and backward directions, as 
presented in the Fig. 16.101 (these are results of calculations performed with 
30 times longer lifetime of Delta resonance and the cross section for its ab- 
sorption enlarged by a factor 1.5). The improvement of agreement between 
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Table 6.2: Multiplicity of pions from p+Au reaction at 0.73 GeV and 2.5 
GeV proton beam energy, without any corrections and with enlarged cross 
section for the absorption of Delta resonance in reaction NA — > NN and 
enlarged Delta resonance lifetime 





Tp = 0.73 GeV 


Tp = 2.5GeV 


without corrections 


with corrections 


without corrections 


with corrections 


7r+: 


0.15 


0.073 


0.46 


0.43 




0.11 


0.053 


0.54 


0.48 


7r~: 


0.065 


0.024 


0.49 


0.43 



the data and model results is substantial, however, the introduced modifi- 
cations are rather drastic. (Improvement of pion dynamics by modifications 
of Delta resonance lifetime and cross section for Delta absorption has been 
also performed in frame of the INCL model in [58j , but obtained results also 
do not comply well with experimental data). Furthermore, it is tested, that 
the influence of the joined corrections on the pion spectra depends on the 
value of proton beam energy. The higher beam energy, the lower reduction 
of the pions yield by the corrections, as it is seen in Figs. 16.111 and I6.12[ for 
spectra of pions produced in p+Au reaction, at low (0.73 GeV) and higher 
(2.5 GeV) beam energy, respectively. The quantitative distinction of the pi- 
ons yield from the reactions at two different incident energies are collected 
in Table 16.21 It has been also verified that implemented corrections do not 
modify yield and shape of spectra of nucleons. Moreover, the corrections do 
not influence the properties of residual nuclei after the first stage of reaction. 
In particular, they do not change excitation energy of the remnants, because 
the average energy carried out by emitted pions with and without the cor- 
rections is the same; due to the corrections less pions are emitted, but they 
carried, in average, more kinetic energy. 

The not very satisfactory agreement between calculations and experimen- 
tal data is not too surprising, because pion production in proton - nucleus 
collisions is a complicated process. Although the elementary cross sections 
for pions production due to free nucleon - nucleon interactions are measured 
and very well known (see e.g. [59]), in-medium problem appears. In the HSD 
model, the in-medium effects are partly included through e.g. employment of 
the Pauli principle and cross section for Delta absorption. However, proper 
handling of modifications due to in-medium effects are hampered by a lack 
of experimental and theoretical knowledge. 
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Pion production mechanism is highly sensitive to pion and A - resonance 
dynamics. So, apart from tested above modifications of the A parameters, 
also introduction of density dependent width of A - resonance (increasing 
with density), introduction of pion potential, change of an angular distribu- 
tion in order to describe possible anisotropic A decay into vrN channel, etc. 
can infiuence the pions production. In fact, till now none of the trials have 
given satisfactory results, see [201 EI]- The propriety of description of pions 
dynamics and production due to proton - nucleus reaction could not be really 
verified, because experimental data are rather scarce, especially in projectile 
energy range of a few GeV order. 

There exists also one more reason that infiuences the calculated pion spec- 
tra. In the HSD model, a target nucleus, i.e. the nuclear density distribution 
is constructed up to a constant limit R + Rcut, where: R = 1.12 ■ y4^/^[/m] 
is a radius of nucleus, defined as a distance at which the value of nucleus 
density decreases to a half of the value of density in the center of nucleus; A 
is a mass number of the nucleus; Rcut = 1.5/m. This is due to difficulties 
with modelling of the infinite tail of nuclear density at nucleus surface, see 
|62j . Although the truncated part corresponds to only a small fraction of the 
whole nucleus density distribution, it can infiuence the pion spectra. This is 
because only pions produced close to the surface have a chance to be emitted, 
as the mean free path of pions is equal to about 1.0 - 1.5 fm, at standard 
density. The introduced cut off of nuclear surface neglects the possibility 
of further pion - nucleon interactions at surface, which are quite probable 
(since pion is a strongly interacting particle) and lead to spuriously higher 
multiplicity of emitted pions. The nucleus density cut off does not influence 
nucleons production, because most of nucleons are emitted closer to the cen- 
ter of nucleus, due to the fact that their mean free path is equal to about 
3.0 fm. Therefore, the HSD model description of proton induced spallation 
reaction is not falsifled by the used limited distribution of density of target 
nuclei, only outgoing pion distributions are to some extent misshapen. 
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Figure 6.3: Two-dimensional distributions of transversal versus longitudinal 
momentum of neutral pions emitted during proton induced reaction on Ni 
target, at various incident energies; results of the HSD model calculations 
(the most central contour line corresponds to a maximum of distribution) 
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Figure 6.5: Multiplicity (left) and average values of kinetic energy per in- 
dividual pion (right ) of pions emitted during proton induced reaction on Al 
target, in function of incident energy; results of the HSD model calculations 
(error bars indicate values of standard deviation of the average values, di- 
vided by a square root of number of events) 
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Figure 6.7: Double differential negative pion spectra from p+Pb reaction at 
0.73 GeV proton beam energy; lines show results calculated with HSD model, 
symbols indicate the experimental data 157] 
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Figure 6.9: Double differential negative pion spectra from p+Pb reaction at 
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Figure 6.10: Double differential negative pion spectra from p+Pb reaction 
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cross section for the absorption of Delta resonance in reaction NA — > NN 
and enlarged Delta resonance lifetime 
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Figure 6.12: Pions kinetic energy spectra from p+Au reaction at 2.5 GeV 
proton beam energy; solid lines show results calculated with HSD model, 
without any corrections, dashed lines show results calculated with enlarged 
cross section for the absorption of Delta resonance in reaction NA — > NN 
and enlarged Delta resonance lifetime 



Chapter 7 



Statistical emission of particles 
in the second stage of proton - 
nucleus collisions 

The second stage of proton - nucleus reaction can be theoretically described 
with various statistical models. Historical origin of the models is related 
to Weisskopf's paper published in 1937 [1]. Weisskopf considers emission 
(evaporation) of a neutron with kinetic energy from nucleus A, excited 
to energy Ea, transforming it to nucleus B, with excitation energy Eb = 
EA — Bn — Ek- Probability of such emission is calculated using the principle of 
detailed balance. It can be determined if cross section for the inverse process 
a{EA, Ek) is known (i.e. cross section for production of compound nucleus 
by collision of a neutron with nucleus B). The probability is expressed as: 

Wn{Ek) = a{EA, Ek)gmEkUJB{EB)/7i'h'uA{EA) (7.1) 

where: m stands for neutron mass, g is a. spin degeneracy factor, uja and ujb 
are the level densities of the initial and final nuclei, respectively. 
In the Weisskopf expression, angular momentum effects are not taken into ac- 
count. In order to consider the effects, one needs to use the Hauser-Feshbach 
formalism ^3]. In this case, the decay probability is determined by level 
density p{Ef, Jj) of the residual nucleus, depending both on its excitation 
energy Ef and angular momentum Jj, and transmission coefficient of emitted 
particle Ti{Ek) {I is the orbital angular momentum removed by the particle, 
Ek is its kinetic energy). 

Decay width F for state of energy E and angular momentum J into a specific 
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channel i is expressed as: 

pE-BEi 

Ti{E, J) = D{E, J)/27T J2 Sil+l+Jf-J) / Ti{Ek)p^iE-BEi-Ek, Jf)dEk 

i,i,Jf ° 

(7.2) 

where: D[E, J) is the level spacing (reciprocal of the level density) of the 
decaying nucleus, / stands for the particle spin. The following energy and 
angular momentum conservation constraints are employed: 
J = I + l + Jf, 

E = Ef + Ek + BEi, where BEi is the binding energy of the emitted particle 
in decay channel i. 

The decay probability into channel i is evaluated as: Ti/Ttot, where: Ttot = 
X^iTj, i runs over all open channels. 




Figure 7.1: The Statistical Evaporation Process 



Important step to comprehensive understanding of the dependence of the 
level density on angular momentum of compound nucleus was done by J. 
R. Grover [M]. He invented the "Yrast level" term to describe the lowest 
lying level at some given angular momentum. It means, that for any angular 
momentum J the nucleus has to have an excitation energy at least equals 
to the rotational energy Erot{J)- Therefore, there cannot be any states with 
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excitation energy below Erot{J), i-e. below "Yrast level". Particle decay 
is inhibited close to the "Yrast line" (i.e. the level density is close to zero 
for the final states allowed by energy and angular momentum conservation). 
Deexcitation proceeds mainly by gamma emission along the "Yrast line", as 
illustrated in Fig. 17.11 

In the frame of this work, calculations of the second stage of proton - 
nucleus reactions have been done with use of PACE2 [20] and GEM 
codes. Below, basic features and results of the models are presented. 

7.1 Evaporation model - PACE 2 

The statistical model code PACE2 (Projection Angular- Momentum Coupled 
Evaporation, version 2) [20] is based on the Hauser-Feshbach formalism. The 
decay sequence of an excited nucleus is calculated with use of a Monte Carlo 
procedure, which is followed until the nucleus reaches its ground state. At 
each stage of deexcitation, the decay probability is determined by calculating 
decay width for neutron, proton, alpha and gamma. Final states differ in 
excitation energy, angular momentum, nucleon number A and proton number 
Z. 

The model contains several parameters, which has to be fixed. The main 
ingredients of the code are the level density and transmission coefficients. 
The level density p{E, J), at a specific value J, is given by: 

p{E, J) = (2 J + l)uj{E - ErotiJ)) (7.3) 

where: uj{E) is the total level density at excitation energy E, evaluated 
according the following formula (given by Weisskopf for Fermi g Ref. 

i, my- 

uj{E) ^ exp[2{aE'f^] (7.4) 

with level density parameter a = const (its dimension is MeV~^). 

For selecting the rotational energy contribution, Erot{J), model of A. J. Sierk 

[UU] is used. 

The transmission coefficients T; are calculated with use of default potentials; 
for protons and neutrons taken from Percy and Percy [6^, for alpha particles 
from Igo and Huizenga [68]. As the nucleus decays, it is assumed that trans- 
mission coefficients of charged particles are shifted in their kinetic energy 
dependence: Ej^ — Vc{A, Z), where A and Z are the nucleon and proton num- 
bers of the emitting nucleus, and Vc is the Coulomb barrier for the emitted 
particle. 
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7.2 Generalized Evaporation Model 



Table 7.1: Multiplicity of particles evaporated from excited Au nucleus, with 
excitation energy equal to 100 MeV and different values of angular momen- 
tum; results of tie PACE2 model calculations 





h 


10 h 


20 h 


n: 


7.80 


7.80 


7.70 


p: 


0.076 


0.080 


0.070 


^He: 


0.031 


0.031 


0.031 



Fission of the excited nuclei is also taken into account. The fission barrier of 
Sierk is used [66] . 

Gamma emission is considered statistically using El, E2, Ml and M2 
transition intensities. 

7.2 Generalized Evaporation Model 

Another model used for the Monte Carlo simulations of the second stage of 
the spallation reaction is the Generalized Evaporation Model (GEM), devel- 
oped by S. Furihata [21], based on the Weisskopf-Ewing approach [UEH]. In 
the GEM model, 66 nuclides up to Mg are included as ejectiles, not only the 
dominant particles emitted from an excited nucleus (i.e. nucleons and helium 
isotopes), as it is in the PACE2. The accurate level density function is used 
for deriving the decay width of particle emission, instead of an approximate 
form, as used in the PACE2 model (see the previous Section). Additionally, 
depending on the excitation energy and mass of the decaying nucleus, differ- 
ent parametrizations of the level density p{E) are applied, see Ref. [2T] . 
In the GEM model, contrary to the PACE2, the dependence of the level 
density on the angular momentum is neglected. Nevertheless, negligible dif- 
ference between spectra of evaporated particles, calculated with the PACE2 
model, with zero and non-zero values of angular momentum of emitting nu- 
clei, has been noticed, as displayed in Fig. 17.21 Adequate multiplicities of 
evaporated particles, calculated with the PACE2 model, with different values 
of angular momentum of emitting nucleus are collected in Table 17. 1[ 
Spectra of evaporated particles (nucleons and ''He) calculated with the PACE2 
and GEM models are in good agreement, as shown in Fig. 17.31 for an exam- 
ple Au nucleus, with the excitation energy set to 100 MeV and zero value 
of angular momentum. Adequate multiplicities of the emitted particles are 
collected in Table [721 
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Table 7.2: Multiplicity of particles evaporated from excited Au nucleus, with 
excitation energy equal to 100 MeV; results of the GEM and the PACE2 
model calculations 





GEM 


PACE2 


n: 


8.86 


7.80 


P: 


0.17 


0.076 


^He: 


0.029 


0.031 



Results of the GEM model calculations are very sensitive on the used level 
density parameter a [21]. There is a possibility to use a parameter equal to 
a quotient of mass of emitting nucleus and a constant number a = A/ const 
MeV~^ (see formula 17.41) . where const ^ 5 ^ 20 (see Ref. [HI EOj) or one 
can use the precise Gilbert - Cameron - Cook - Ignatyuk (GCCI) level den- 
sity parametrization with the energy dependence taken into account, see 
Ref . [2T1 171] . Dependence on the used method of level density parametriza- 
tion is noticeable e.g. in multiplicities and spectra of evaporated particles. 
Comparison of spectra of various isotopes evaporated during the second stage 
of an example p+Au reaction, at 2.5 GeV beam energy, results of the GEM 
model calculations, with the different level density parameters are presented 
in Figures 17.41 and 17.51 It is seen from the Figures, that distributions of 
light particles (isotopes of hydrogen and helium) calculated with different 
level density parameters are very similar. Contrary, distributions of heav- 
ier isotopes (lithium, beryllium, etc.) differ significantly. Use of the GCCI 
level density parameter gives much higher abundance of the heavier isotopes, 
than use of the simple a = A/8 parameter. Moreover, calculations with the 
a = A/20 parameter result with even higher abundance than with the GCCI 
level density parameter. Adequate multiplicities of the evaporated particles 
are collected in Table 17. 3[ It is seen from the Figures, that choice of level 
density parameter affects not only heights of the spectra, but influences also 
their slopes. In the frame of the work, the GCCI level density parameter is 
used (confrontation with experimental data leads to conclusion, that use of 
the GCCI level density parameter gives the most reliable results). 
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7.2 Generalized Evaporation Model 



Table 7.3: Multiplicity of particles evaporated during p+Au reaction, at 2.5 
GcV proton beam energy, calculated with the GCCI, a — A/S and a — A/20 
level density parameter 





GCCI 


a^A/S 


a = A/20 


n: 


14.78 


14.90 


9.016 


p: 


1.76 


1.97 


0.85 


^H: 


1.11 


1.11 


1.17 


3H: 


0.51 


0.44 


0.78 




0.063 


0.062 


0.12 


7 

^He: 


0.84 


0.76 


0.85 


^He: 


0.012 


0.0083 


0.059 


^Li: 


0.020 


0.013 


0.17 


— ^tT 

^Li: 


0.018 


0.012 


0.22 


«Li: 


0.0055 


0.0032 


0.104 


9Li: 


0.0011 


0.00057 


0.035 


^Be: 


0.0021 


0.0016 


0.027 


9Be: 


0.0039 


0.0024 


0.081 


lOBe: 


0.0045 


0.0024 


0.11 


iiBe: 


0.00041 


0.00019 


0.021 


8B: 


0.00006 


0.00004 


0.0013 


lOB: 


0.0022 


0.0013 


0.060 


"B: 


0.0018 


0.00091 


0.076 


12B: 


0.0011 


0.00053 


0.039 


13B: 


0.00013 


0.00005 


0.0049 


C: 


0.0027 


0.0015 


0.093 


N: 


0.0011 


0.00055 


0.071 


0: 


0.00041 


0.00025 


0.034 
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Figure 7.2: Spectra of particles evaporated from excited Au nucleus, with 
the excitation energy equal to 1 00 MeV and different values of angular mo- 
mentum; results of the PACE2 model (solid lines - zero value of angular 
momentum of the excited nucleus, dashed lines - angular momentum equal 
to lOh, dotted hnes - 20 h) 
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Figure 7.3: Spectra of particles evaporated from excited Au nucleus, with the 
excitation energy equal to 100 MeV and zero value of angular momentum; 
solid lines - results of the PACE2 model, dashed lines - results of the GEM 
model 
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Figure 7.4: Spectra of isotopes evaporated during p+Au reaction, at 2.5 GeV 
beam energy; results of the HSD+GEM model calculations (solid lines - cal- 
culations with the GCCI level density parameter, dashed lines - calculations 
with a^A/8, dotted lines - a = A/20) 
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7.2 Generalized Evaporation Model 
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Figure 7.5: Spectra of isotopes evaporated during p-hAu reaction, at 2.5 GeV 
beam energy; results of the HSD+GEM model calculations (solid lines - cal- 
culations with the GGGI level density parameter, dashed lines - calculations 
with A/8, dotted lines - a = A/ 20) 
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Bulk models predictions for the 
proton induced nuclear 
reactions 



8.1 Illustration of energy balance of the re- 
action 

In order to show, how large fraction of energy is carried out during the first 
and the second stage of proton - nucleus reaction, respectively, a sample case 
of the HSD and statistical evaporation model calculations, for p + Au colli- 
sion at 2.5 GeV proton beam energy is analysed in this Section. 
The total initial kinetic energy introduced into nucleus by incoming proton 
is equal to 2.5 GeV. During the first stage of the reaction, about 1.7 GeV of 
the energy is carried away by nucleons and about 0.4 GeV - by pions. The 
multiplicity of nucleons is near 4.43, whereas the pion multiplicity, summed 
over all pion charges, is equal to 1.48. The average kinetic energy of one 
nucleon is about 0.4 GeV, whereas the average kinetic energy of one pion 
is equal to 0.29 GeV. Non-negligible amount of energy is taken away from 
the excitation energy and converted into rest masses of emitted pions, i.e. 
0.2 GeV. Remaining energy of the excited nucleus after the first stage of the 
reaction, is equal to only about 200 MeV. This amount of energy is then 
evaporated in the second stage of the reaction by emission of mainly neu- 
trons, protons, alphas and gammas. Emitted neutrons accumulate about 20 
MeV of the energy, protons carried about 3 MeV, alphas - about 1 MeV, 
and gammas - about 6 MeV. The multiplicity of neutrons is equal to about 
10 and the average kinetic energy of one neutron is equal to about 2 MeV, 
multiplicity of protons is equal to about 1 and the average kinetic energy of 
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8.2 Participation of Gssion processes 



one proton is equal to about 3 MeV. Multiplicity of alphas is equal to about 
0.3 and the average kinetic energy of one alpha is equal to about 4 MeV and 
finally multiplicity of gammas is equal to about 6 and the average kinetic 
energy of one gamma is equal to 1 MeV. After taking into account the bind- 
ing energy of nucleons inside nucleus, i.e. about 8 MeV per one nucleon, the 
whole amount of excitation energy of the residual nucleus (about 200 MeV) 
is carried out during the second stage of the reaction. 

The presented balance shows, that the biggest amount of initial energy de- 
posited in nucleus by incoming proton is carried out during the first stage of 
the reaction, i.e. about 85 %, and only about 15 % of the initial energy is 
carried away during the second stage of the reaction. 

8.2 Participation of fission processes in spal- 
lation reaction 

During the second stage of the proton - nucleus spallation reaction, deexcita- 
tion of the hot residual nucleus by evaporation of various isotopes takes place, 
as it is mentioned above in this work. However, there is a possibility, that 
during this process the evaporating nucleus undergoes fission. In such case, 
the particles emission is continued from the fission fragments. Probability 
of fission of particular nuclei depends on the ratio of Coulomb energy and 
surface energy of individual nucleus. The probability is strongly dependent 
on fissility parameter (/) defined as: 



where: Z is the charge of a nucleus and A is the mass number of a nucleus. 

For example, the fissility parameter for Au nucleus, calculated according 
to the definition (18.11) is equal to 31.7. The experimental data collected in 
Ref. [72] indicate, that Au nucleus has rather low probability for fission, in 
comparison of Th and U nuclei. Distribution of the fissility parameter for 
residual nuclei emerged after the first stage of an example p+Au reaction, at 
2.5 GeV beam energy (result of the HSD model calculations) is presented in 
Fig. 18. 1[ The average value of the fissility parameter distribution is equal to 
31.40 ±0.97. 

On the base of a behavior of the fissility parameter in function of mass 
of residual nuclei, it can be tested if a ratio of charge and mass number is 
constant during the reaction and equal to the ratio of charge and mass num- 
ber of the initial target nucleus. The fissility parameter in function of mass 
of residual nuclei for the example p+Au reaction, at 2.5 GeV proton beam 
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Figure 8.1: Distribution of the hssility parameter for residual nuclei emerged 
after the first stage of p+Au reaction, at 2.5 GeV beam energy; results of 
the HSD model calculations 
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Figure 8.2: The fissility parameter in function of mass of residual nuclei for 
the example p+Au reaction, at 2.5 GeV beam energy; results of the HSD 
model calculations (x^ means chisquare per degree of freedom ) 



energy is presented in Fig. 18.21 The ratio of charge and mass number of the 
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Au initial target nucleus is equal to: Z(Au)/A(Au) = 79/197 = 0.40. Using 
the definition of the fissility parameter fl8.ll) and assuming that: Z=0.40-A, 
it is obtained that: f=0.16-A. Parametrization of the dependence presented 
in the Fig. 18.21 by a function proportional to the mass number give exactly 
the same result as received from the definition. This indicates, that in a 
first approximation, a ratio of charge and mass number is constant during 
the reaction and equal to the ratio of charge and mass number of the initial 
target nucleus. However, parametrization of the dependence by a function 
proportional to some power of the mass number (correcting a slope) gives 
a bit better accuracy, as presented in Fig. 18.31 It means, that the ratio of 
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Figure 8.3: The fissility parameter in function of mass of residual nuclei for 
the example p+Au reaction, at 2.5 GeV beam energy; results of the HSD 
model calculations (x^ means chisquare per degree of freedom) 
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charge and mass number is not exactly equal to the ratio for the initial target 
nucleus. This is confirmed by multiplicities of protons and neutrons emitted 
during first stage of the p+Au reaction, at 2.5 GeV proton beam energy. 
Namely, due to interactions inside the target nucleus, multiplicity of emitted 
protons (equal to 1.66) is a bit higher than multiplicity of emitted neutrons 
(equal to 1.63). Contrary, in case of n+Au reaction, at 2.5 GeV neutron 
beam energy, multiplicity of emitted neutrons during the first stage of the 
reaction (equal to 1.89) is a bit higher than multiplicity of emitted protons 
(equal to 1.38). This indicates, that information about type of projectile is 
kept during the first stage of reaction. 

If looking at the kinetic energy distributions of isotopes emitted during 
the second stage of p+Au reaction, calculated with blocked and allowed fis- 
sion, differences are visible, because Au nucleus possesses non-zero, although 
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rather small fissility. The example spectra (result of the GEM model calcu- 
lations) are presented in Fig. 18. 4[ It is clearly seen from the Figure, that 
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Figure 8.4: The kinetic energy distributions of ^He evaporated during the 
second stage of p+Au reaction, at 2.5 GeV beam energy; results of the 
HSD+GEM model calculations (solid line - fission allowed, dashed line - 
fission blocked) 



fission shifts spectra of emitted particles down in their maximum and causes 
them wider in energetical range, in direction of both lower and higher ener- 
gies (in Fig. 18.41 the effect is shown for alpha particles). This is because of 
evaporation from fission fragments, which move in opposite directions. As a 
consequence there is a bit less particles with intermediate energies, emitted 
from the unfissioned residuum, but more particles with a bit lower and higher 
energies, emitted after fission, from excited fragments moving in backward 
and forward directions, respectively. 

Cross section for the different fragments production in function of the 
fragments masses and charges (result of the HSD+GEM model calculations), 
for the example p-l-Au reaction, at 2.5 GeV proton beam energy, calculated 
with blocked and allowed fission, are presented in Fig. 18.51 Interesting result 
is obtained, when one subtracts the cross sections obtained with blocked and 
allowed fission. One can see in the Figures that fission diminishes amount of 
nuclei of masses close to 150 - 160 and increases amount of nuclei of masses in 
the range 60 - 100. One can conclude, that fission is much slower process than 
evaporation. First evaporation from a hot residual nucleus takes place. Due 
to evaporation both the mass and excitation energy of residuum decrease. 
Also the Coulomb energy and the surface energy of the nucleus is changing. 
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If the emission of particles is slow, a possibility of fission occurs. Fission is 
a collective slow process. This is confirmed by experimental data collected 
in Ref. [72j, which show, that fission cross section decreases with increasing 
energy provided into nuclei. In case of very fast evaporation, fission has no 
possibility to occur. The second stage of spallation reaction is a competition 
between evaporation and fission. 
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Figure 8.5: Cross section for different fragments production in function of 
fragment mass and charge from p+Au reaction, at 2.5 GeV beam energy, 
with (a) allowed and (b) blocked fission; results of the HSD+GEM model 
calculations; (c) differences between calculations with and without fission 
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Comparison of results of the 
HSD plus evaporation model 
calculations with experimental 
data 

Result of the HSD plus statistical evaporation model calculations have been 
examined by confrontation with existing experimental data. Effects are dis- 
cussed in this Chapter. 

9.1 Neutron spectra 

Inclusive differential kinetic energy spectra of neutrons, emitted at different 
angles during proton - nucleus reactions on various target nuclei and at var- 
ious impact energies, calculated with the HSD plus evaporation model and 
compared with available experimental data are analysed in this Section. 
High energy part of neutron distribution, composed of so-called fast neu- 
trons, produced during the first stage of the reaction is a result of the HSD 
model calculations. Such neutrons are emitted mainly in forward direction. 
Low energy part of neutron distribution, composed of slow neutrons produced 
during the second stage of the reaction is a result of evaporation model calcu- 
lations. Such neutrons are emitted both in forward and backward directions. 
Neutron spectra obtained as a result of summing up these two parts are 
compared with the following experimental data. Distributions of spallation 
neutrons produced in proton induced reactions on Al. Fc. Zr, W, Pb and Th 
targets at 1.2 GcV beam energy and on Fe and Pb targets at 1.6 GeV proton 
beam energy, emitted at angles: 10°, 25°, 40°, 55°, 85°, 100°, 115°, 130°, 145° 
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and 160° are compared with data measured at the SATURNE accelerator in 
Saclay [73]. Similarly, calculated spectra of neutrons produced in proton in- 
duced reactions on Pb target at 0.8 and 1.5 GeV beam energies, emitted at 
angles: 15°, 30°, 60°, 90°, 120° and 150° are compared with experimental 
data measured by Ishibashi et al. \T^. The comparisons are presented in 
Figures: ElEIllin^linainSlinSlEIliniHliniiand ICTl respectively. 

One can see, that independently on the value of projectile energy, the 
calculated neutron spectra produced in proton induced reactions on light 
targets (i.e. Al, Fe and Zr, presented in the Figures 19.11 19.21 19.31 and 19.71 
respectively) are in perfect agreement with the adequate experimental data. 
Contrary, in case of spectra of neutrons produced in proton induced reac- 
tions on heavy targets (i.e. W, Pb and Th, presented in the Figures 19.51 
19. 6119. 8t 19.91 and [5. lot respectively), the good agreement between calculations 
and proper experimental data is seen only in backward angles. In forward 
angles, calculations underestimate experimental data in the central parts of 
the distributions. Furthermore, based on the presented results of the several 
reactions, on the six target nuclei, the regular growth of the missing part 
of the calculated spectra is observed. The effect is continuous and increase 
monotonically with mass of target. The missing parts, in case of reactions 
on heavy targets indicate, that apart from the two stages of spallation reac- 
tion (fast and slow stage), there must be an additional intermediate stage, 
probably preequilibrium stage. It seems, that neutron production during this 
intermediate stage is dominant in case of reactions with heavy targets and 
negligible for reactions with light targets. 

Multiplicities of neutrons emitted during first and second stage of the 
reactions are collected in Table 19. 1[ It is clearly seen, that neutrons are 
emitted mainly in the second stage of spallation reaction. 
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Table 9.1: Multiplicities of neutrons emitted during fast and slow stage of 
reaction 
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Figure 9.1: Double differential neutron production cross section for p+Al re- 
action, at 1.2 GeV proton beam energy; lines show results of the HSD+GEM 
model calculations (dashed and dotted lines are contributions of first and 
second stage of reaction, solid line is their sum), symbols indicate the exper- 
imental data 
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Figure 9.2: Double differential neutron production cross section for p+Fe re- 
action, at 1.2 GeV proton beam energy; lines sliow results of the HSD-hGEM 
model calculations (daslied and dotted lines are contributions of first and 
second stage of reaction, solid line is their sum), symbols indicate the exper- 
imental data 1731 



Chapter 9. Comparison of calculations with experimental data 111 



1 oo 
1 o 
1 

0.1 
0.01 
0.001 







1 o° 











0.001 0.01 0.1 1 

Er, (OeV) 






0.001 0.01 



n-nX T„=1 .2 GeV. 





0.001 0.01 




0.001 0-01 




0.001 0.01 




p-^ ^^Zr n+X Tp=1 .2 GeV 




Figure 9.3: Double differential neutron production cross section for p+Zr re- 
action, at 1.2 GeV proton beam energy; lines show results of the HSD+GEM 
model calculations (dashed and dotted lines are contributions of first and 
second stage of reaction, solid line is their sum), symbols indicate the exper- 
imental data fTSl 
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Figure 9.4: Double differential neutron production cross section for p-hW re- 
action, at 1.2 GeV proton beam energy; lines sliow results of the HSD-hGEM 
model calculations (daslied and dotted lines are contributions of first and 
second stage of reaction, solid line is their sum), symbols indicate the exper- 
imental data 1731 
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Figure 9.5: Double differential neutron production cross section for p-hPb re- 
action, at 1.2 GeV proton beam energy; lines show results of the HSD+GEM 
model calculations (dashed and dotted lines are contributions of first and 
second stage of reaction, solid line is their sum), symbols indicate the exper- 
imental data fTSl 
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Figure 9.6: Double differential neutron production cross section for p+Th re- 
action, at 1.2 GeV proton beam energy; lines show results of the HSD+GEM 
model calculations (dashed and dotted lines are contributions of first and 
second stage of reaction, solid line is their sum), symbols indicate the exper- 
imental data 1731 



Chapter 9. Comparison of calculations with experimental data 115 




o 




0.001 

0.001 




0.001 

0.001 




o 



En (GeV) En (GeV) 




O 



En (GeV) En (GeV) 




0.001 

0.001 



1 O 



1 OO 
1 O 
1 

O.I 
0.01 



0.001 

0.001 



p-i- Fe — 


^ n-i-X Tp=1 .6 GeV 








1 60° 







0.01 O.I 1 
En (GeV) 



1 O 



Figure 9.7: Double differential neutron production cross section for p+Fe re- 
action, at 1.6 GeV proton beam energy; lines show results of the HSD+GEM 
model calculations (dashed and dotted lines are contributions of first and 
second stage of reaction, solid line is their sum), symbols indicate the exper- 
imental data fTSl 



llF 



116 

1 ooo 
1 oo 
1 o 
1 

0.1 
0.01 
0-001 





->■ n-i-X X_ = 1 


e QeV 






M 

1 o° 



















0.001 0.01 



0.1 

(aeV) 



9.1 Neutron spectra 



1 ooo 
1 oo 
1 o 
1 

0.1 
0.01 
0.001 





^ rn-X Xp=1 


.6 QeV 




25° 


\ 
1 


Y 





0.001 0.01 



0.1 1 
(aeV) 



1. 




0.001 0.01 




0.001 0.01 




0.001 

0.001 0.01 






0.001 0.01 




0.001 0.01 




0.001 0.01 



Figure 9.8: Double differential neutron production cross section for p-hPb re- 
action, at 1.6 GeV proton beam energy; lines show results of the HSD+GEM 
model calculations (dashed and dotted lines are contributions of first and 
second stage of reaction, solid line is their sum), symbols indicate the exper- 
imental data 1731 
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Figure 9.9: Double differential neutron production cross section for p-hPb re- 
action, at 0.8 GeV proton beam energy; lines show results of the HSD+GEM 
model calculations (dashed and dotted lines are contributions of first and 
second stage of reaction, solid line is their sum), symbols indicate the exper- 
imental data f7^ 
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Figure 9.10: Double differential neutron production cross section for p+Pb 
reaction, at 1.5 GeV proton beam energy; lines show results of the 
HSD-hGEM model calculations (dashed and dotted lines are contributions of 
first and second stage of reaction, solid line is their sum), symbols indicate 
the experimental data \74]j 
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9.2 Proton spectra 

Inclusive, double differential kinetic energy spectra of protons emitted at 
different angles in proton induced reactions on various targets, at broad range 
of incident energy, calculated with the HSD+GEM model, are compared 
with available experimental data. Model results of low energy reactions: 
p+Ni at 0.175 GeV beam energy (in this case, time of first stage calculations 
is equal to 55 fm/c), compared with experimental data measured by PISA 
collaboration [7j and p+Bi at 0.45 GeV beam energy (with time of first 
stage calculations equals to 65 fm/c), compared with data published in [73], 
are presented in Figures 19.111 and 19.121 respectively. Very good agreement 
between calculations and data is clearly seen. 

Comparisons of calculated proton spectra produced in reactions induced 
by higher incident energies, e.g. p+Au at 1.2 GeV, 1.9 GeV and 2.5 GeV 
projectile energy (with time of first stage calculations equal to 45 fm/c), con- 
fronted with data measured by PISA collaboration [7], are shown in Figures 
19.131 19.141 and 19.151 In order to make the comparison easier, the high en- 
ergy parts of the calculated distributions are plotted only up to 200 MeV of 
proton kinetic energy. It is seen, that in case of 1.2 GeV of incident energy 
(Fig. I9.13p . the high energy parts of experimental distributions are described 
well, the low energy parts are slightly overestimated by calculations. In the 
case of 1.9 GeV beam energy (Fig. 19.141) . the whole spectra, especially pro- 
tons emitted in backward directions are described very well. Low energy 
part of experimental distributions of protons from reaction at 2.5 GeV (Fig. 
I9.15P are also very well described by calculations. The high energy parts are 
slightly underestimated. Nevertheless, the average agreement between the 
HSD+GEM model results concerning proton spectra and adequate experi- 
mental data is satisfactory. 

In the presented Figures, contributions from intra-nuclear cascade (high 
energy part, result of the HSD model calculations) and evaporation (low 
energy part, result of the GEM model calculations) are visible. The cascade 
contribution evidently dominates. This indicates, that protons are emitted 
mainly in the first stage of spallation reaction. This is additionally confirmed 
by values of multiplicities of protons emitted during first and second stage of 
reaction, compared in Table 19.21 (quite contrary to the neutron case, consider 
Table EH). 

Nucleon distributions depend strongly on projectile energy. The example 
two-dimensional distributions of transversal (P^) versus longitudinal (Pz) 
momentum of nucleons emitted during the first stage of proton induced re- 
actions on Ni target, at various incident energies are presented in Fig. 19.161 
These are results of HSD model calculations, where Pz corresponds to the nu- 
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Table 9.2: Multiplicities of protons emitted during first and second stage of 
reaction 



reaction 


first stage 


second stage 


p(0.175 GeV)+Ni 


1.61 


0.91 


p(0.45 GeV)+Bi 


2.16 


1.00 


p(1.2 GeV)+Au 


2.67 


1.13 


p(1.9 GeV)+Au 


2.42 


1.099 


p(2.5 GeV)+Au 


2.20 


0.96 



cleons momentum in beam direction and Pt is defined as Pt = y Px + ^y- 
The two-dimensional distributions depict, how the nucleon spectra evolve 
with incident energy. In all the cases, low energy part is symmetric around 
zero value (it is not so in the case of pion production, see Fig. 16. 3p . With 
increase of projectile energy, the distribution is building up mainly in longi- 
tudinal (forward) direction. It spreads always up to a value corresponding 
to the value of beam momentum. That is because, it is possible that incom- 
ing proton will scatter only slightly on a surface of target nucleus and will 
appear in forward direction with almost unchanged momentum. In the Fig. 
19.161 momentum spectra of all nucleons are presented, since distributions 
of protons and neutrons emitted during first stage of the reaction are very 
similar. The only exception concerns momentum corresponding to the beam 
momentum. 
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Figure 9.11: Double differential proton spectra from p+Ni reaction at 0.175 
GeV incident energy; lines show results of HSD and GEM model calculations 
(solid line - first stage, dashed line - second stage), symbols indicate the 
experimental data measured by experiment PISA I7j 
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Figure 9.12: Double differential proton spectra from p+Bi reaction at 0.45 
GeV proton beam energy; lines show results of HSD and GEM model calcu- 
lations (solid line - first stage, daslied line - second stage), symbols indicate 
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Figure 9.13: Double differential proton spectra from p+Au reaction at 1.2 
GeV proton beam energy; lines show results of HSD and GEM model calcu- 
lations (solid line - first stage, dashed line - second stage), symbols indicate 
data measured by experiment PISA 
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Figure 9.14: Double differential proton spectra from p+Au reaction at 1.9 
GeV proton beam energy; lines show results of HSD and GEM model calcu- 
lations (solid line - first stage, dashed line - second stage), symbols indicate 
data measured by experiment PISA ITj 
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Figure 9.15: Double differential proton spectra from p+Au reaction at 2.5 
GeV proton beam energy; lines show results of HSD and GEM model calcu- 
lations (solid line - first stage, dashed line - second stage), symbols indicate 
data measured by experiment PISA 
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Figure 9.16: Two-dimensional distributions of transversal versus longitudi- 
nal momentum of nucleons emitted during the first stage of proton induced 
reaction on Ni target, at various incident energies; results of the HSD model 
calculations (the most central contour line, symmetrical around zero value, 
corresponds to a maximum of distribution) 
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9.3 Deuteron spectra 

Comparison of inclusive differential energy spectra of deuterons calculated 
with the HSD+GEM model, with example experimental data is presented in 
Fig. EITl 
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Figure 9.17: Double differential deuteron spectra from p+Au reaction, at 
2.5 GeV proton beam energy; lines show results of the HSD+GEM model 
calculations, symbols indicate data measured by PISA experiment '0 

One can see, that only low energy part (i.e. deuterons evaporated from 
the residual nucleus during second stage of spallation reaction) is described 
by the model calculations. The absence of high energy part of modelled 
distributions indicates a lack of deuterons produced during first stage of the 
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reaction. Certainly, description of the whole distributions necessitates imple- 
mentation of some additional mechanisms into the first stage model. Since 
deuteron spectra behave very similar as proton distributions (compare e.g. 
Figures [9. 151 and [9.17p . one can suppose that the high energy part of deuteron 
spectra are produced due to a coalescence of proton and neutron during first 
stage of reaction. Production of clusters (i.e. light ions) cannot be followed 
with the HSD approach, since it is a one body theory. Nevertheless, a trial of 
construction of deuterons from the fast stage of the spallation reaction, with 
use of informations provided by the HSD approach, have been undertaken 
in the frame of this work. The idea of generation of such deuterons, called 
here deuterons of coalescence, is based on the information of positions and 
momenta of nucleons belonging to the residual nucleus and nucleons being 
emitted during the first stage. The information is taken at a certain point of 
time, i.e. in the end of the first stage calculations. Deuterons can be built 
up in two manners: by fast proton and proper fast neutron, which are al- 
ready emitted, or by fast proton emitted during fast stage of the reaction and 
proper neutron which belongs actually to the excited residual nucleus, but 
could have been glued with and taken away by the escaping proton during 
the first stage of reaction. In both cases, it is assumed that probability of 
creation a deuteron by defined above nucleon pair, depends on their relative 
momentum (Ap). It means, two nucleons are considering as a deuteron if 
their relative momentum fulfills assumed condition. The following conditions 
have been tested: each pair of neutron and proton can create a deuteron with 
the same probability (this gives the highest border of the absolute normal- 
ization), next, deuterons can be created by proton and neutron only if their 
relative momentum is lower than or equal to, respectively: Ap < 0.45 GeV/c, 
Ap < 0.30 GeV/c and Ap < 0.15 GeV/c. 

All formed deuterons are additionally repulsed by Coulomb force. The Coulomb 
potential is approximated by potential of homogeneously charged sphere. 
This is because an incident proton does not cause any significant changes in 
the nucleon density distribution inside target nucleus, as shown in Chapter 
m Therefore, the kinetic energy of each deuteron is enlarged by value: 

Ecoui = ^ [GeV\ (9.1) 

where: Zt and are charges of target nucleus and deuteron, respectively, 
Rt is a radius of target. 

Obtained differential energy spectra of deuterons produced in p+Au reaction 
at 2.5 GeV proton beam energy are presented in Figures [9.181 and 19. 19[ The 
calculated distributions are compared with experimental data measured by 
NESSI collaboration 
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In more precise definition, deuteron is consider as a set of nucleons, wliicli 
are sufficiently close to each other in phase space. Therefore, distributions of 
deuterons constructed of pairs of fast nucleons, with condition set not only 
on their relative momentum (Ap), but also their relative position (Ar) have 
been tested, too. The following unit volume of phase space have been se- 
lected: ApAr < 0.387 GeV-fm/c (= the Fermi momentum x 1.4 fm) [77] . 
The result underestimates experimental data, as shown in Fig. 19.181 and 
19.191 That is because of lost of information about the time evolution of rel- 
ative positions of the nucleons during the reaction. The condition of relative 
positions of nucleons, imposed on the situation at the end of first stage of 
reaction rejects nucleon pairs, that it has no record whether in some earlier 
period of time nucleons were close enough to create deuterons. Therefore, it 
determines, at most, the lowest border of coalescence effect. 
There is also possible to define the probability (in fact the density of proba- 
bility) of creation a deuteron by a pair of nucleons depending on their relative 
momentum, as equal to the modulus of the deuteron wave function squared: 
[■^oIp)!^ "P^- The wave function can be parametrized itself as: 

Mp) = E Aexp{-aip^) (9.2) 

i 

Values of the coefficients Ai and a, are taken from |78j. The used parametriza- 
tion of the deuteron wave function is accurate within a few percent, up to 
momenta of 1 GeV/c, see Ref. |78j . 

Obtained in this way spectra of deuterons underestimate the experimental 
data, see Figures 19. 181 and l 9. 191 

The calculated spectra of deuterons constructed by two fast nucleons (Fig. 
19.181 and 19.201) are normalized absolutely. The spectra of deuterons con- 
structed by fast protons and neutrons from the residual nucleus (Fig. I9.19P 
cannot be normalized absolutely. It is because the information, whether a 
fast proton before leaving target nucleus has been near a neutron long enough 
to create a deuteron is lost, what gives in this case evident overestimation of 
experimental data. 

It is seen from the Fig. 19.201 that the high energy part of the experi- 
mental spectra are well described by deuterons of coalescence created by a 
pair of fast nucleons with relative momentum lower or equal to about 0.40 
GeV/c. It is also seen that distributions of deuterons of coalescence are per- 
fect continuation of the low energy part. As result, it gives full description 
of experimental data. However, taking into account shapes of the distribu- 
tions from Fig. 19. 19^ one cannot reject possibility that some of the deuterons 
are produced as coalescence of e.g. fast proton and neutron from residual 
nucleus. 
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Figure 9.18: Differential spectra of deuterons of coalescence of proton and 
neutron emitted during fast stage ofp+Au reaction, at 2.5 GeV proton beam 
energy, compared with experimental data measured by NESSI j76| ("upper 
limit" means all proton - neutron pairs convert into deuterons) 
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Figure 9.19: Differential spectra of deuterons of coalescence of fast emitted 
proton and a neutron from residual nucleus remained after fast stage ofp+Au 
reaction, at 2.5 GeV proton beam energy, compared to NESSI data [76^ 
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Figure 9.20: Double differential deuteron spectra from p+Au reaction, at 
2.5 GeV proton beam energy; lines sliow results of tie HSD+GEM model 
calculations, symbols indicate the experimental data measured by NESSI \76l 
("GEM" means "statistical evaporation only") 
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9.4 Other Ejectiles 

Looking at comparison of distributions of heavier particles emitted in proton 
- nucleus reaction, calculated with the HSD+GEM model with experimental 
data, displayed in Fig. I9.2H it is seen that only low energy part is described. 
That is because, in the frame of the model, the light composites are produced 
only in evaporation stage, as discussed in Sec. 19.31 

It is interesting to know, what is the predictive power of the model, how 
distributions of the evaporated composites vary with projectile energy and 
mass of target. Example energy spectra of different ejectiles (isotopes of H, 
He, Li, Be, B, ...) emitted in proton induced reaction on Au target, at a 
few chosen beam energies are presented in Figures 19.221 and I9.23[ Spectra 
of particles ejected during 2.5 GeV proton induced reactions on few example 
targets are shown in Figures 19.241 and 19.251 Looking at the distributions 
of each ejectile, it is clearly seen, that shape of the particular distributions 
almost does not change with projectile energy, but strongly depends on used 
target. In all presented cases, the heavier target, the broader the distribu- 
tions. Average kinetic energy of emitted particle increases with increase of 
mass of used target, but in this case neutrons are an exception. This indi- 
cates, that reason for such behavior is connected with charge of the particles, 
i.e. Coulomb repulsion force between the ejectiles and an emitting source. 
The higher charge of the emitting source, the higher average kinetic energy 
gained by charged ejectiles due to the Coulomb repulsion. Similarly, the 
particles yield almost does not change with incident energy, but varies sig- 
nificantly with mass of target. Multiplicities of the particles, produced due 
to discussed here reactions are compared in Tables 19.31 and 19.41 The yield 
of neutrons increases linearly with increase of the target mass. Multiplicity 
of each individual charged ejectile varies differently with the target mass. In 
some cases, it first decreases with the target mass number (for example pre- 
sented in the Table 19.41 the minimal yield corresponds usually to Ni target) 
and then increases with further increase of mass of target (e.g. proton yield 
behaves in opposite way). In many cases, it is also observed roughly linear 
decrease of the multiplicity with increase of the target mass (as e.g. in case of 
^Li). The trends for the particles yield can be to some extent qualitatively 
understood, since it is known, that the particles emission is influenced by 
the Coulomb barrier, which is higher for heavier nuclei and by the binding 
energy, which in average first increases with mass number (A) of target (up 
to A equal to about 56) and then decreases with A. The relative yield of 
different ejectiles emitted in one chosen reaction ( Tables [9731 and is also 
determined by separation energy of individual ejectile, which is not linear 
with the mass number of ejectile. 
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Figure 9.21: Example double differential spectra of various ejectiles emit- 
ted during an example reaction p+Au, at few incident energies; results of 
the HSD+GEM model calculations (solid lines) compared with experimental 
data measured by PISA at COSY ^ (squares) 
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Figure 9.22: Kinetic energy spectra of different ejectiles emitted in proton 
induced reaction on Au target, at several proton beam energies; results of 
the HSD+GEM model calculations 
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Figure 9.23: Kinetic energy spectra of different ejectiles emitted in proton 
induced reaction on Au target, at several proton beam energies; results of 
the HSD+GEM model calculations 
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Figure 9.24: Kinetic energy spectra of different ejectiles emitted in 2.5 GeV 
proton induced reaction on several targets; results of the HSD-hGEM model 
calculations 
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Figure 9.25: Kinetic energy spectra of different ejectiles emitted in 2.5 GeV 
proton induced reaction on several targets; results of the HSD-hGEM model 
calculations 
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Table 9.3: Multiplicities of various ejectiles emitted during second stage of 
proton induced reactions on Au target, at a few impact energies 



iriij)act energy 


1.2 GeV 


1.9 GeV 


2.5 GeV 


n: 


14.17 


12.79 


11.88 


^H: 


1.13 


1.099 


0.96 




0.59 


0.59 


0.51 


^H: 


0.28 


0.28 


0.24 


^He: 


0.020 


0.023 


0.021 


^He: 


0.59 


0.56 


0.48 


^He: 


0.0054 


0.0058 


0.0052 


^Li: 


0.0045 


0.0058 


0.0056 


''Li: 


0.0067 


0.0076 


0.0067 


'Li: 


0.0018 


0.0022 


0.0019 


'Li: 


0.00039 


0.00045 


0.00039 


^Be: 


0.00045 


0.00060 


0.00055 


^Be: 


0.0015 


0.0017 


0.0015 


in T-> 


0.0018 


0.0020 


0.0017 


11 D^. 


0.0001 / 


0.0001 / 


n nnn 1 
0.00016 


'B: 


0.000010 


0.000018 


0.000011 


105: 


0.00065 


0.00081 


0.00069 


"5: 


0.00068 


0.00078 


0.00067 


1^5: 


0.00047 


0.000497 


0.00043 


135: 


0.000061 


0.000051 


0.000048 


C: 


0.0011 


0.0012 


0.0012 


N: 


0.00044 


0.00045 


0.00043 


0: 


0.0002 


0.00021 


0.00017 
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Table 9.4: Multiplicities of various ejectiles emitted during second stage of 
2. 5 GeV proton induced reactions on several targets 



reaction 


P+Al 


P+Ni 


p+Ag 


p+Au 


n: 


0.97 


1.76 


6.11 


11.88 


^H: 


0.82 


2.20 


2.19 


0.96 


^H: 


0.36 


0.51 


0.75 


0.51 




0.076 


0.089 


0.20 


0.24 


^He: 


0.082 


0.10 


0.082 


0.021 


^He: 


0.58 


0.53 


0.64 


0.48 


^He: 


0.0017 


0.0012 


0.0036 


0.0052 


^Li: 


0.033 


0.021 


0.017 


0.0056 


''Li: 


0.0097 


0.0057 


0.0092 


0.0067 


^Lt: 


0.0012 


0.00065 


0.0016 


0.0019 




0.00016 


0.000046 


0.00020 


0.00039 


^Be: 


0.0087 


0.0059 


0.0032 


0.00055 


^Be: 


0.0040 


0.0013 


0.0019 


0.0015 




0.0016 


0.00055 


0.0013 


0.0017 




0.000075 


0.000015 


0.000068 


0.00016 


^B: 


0.00059 


0.00046 


0.00015 


0.000011 




0.013 


0.0029 


0.0019 


0.00069 




0.0061 


0.00088 


0.00079 


0.00067 


i2fi: 


0.00087 


0.00016 


0.00030 


0.00043 


135: 


0.000092 


0.0000069 


0.000018 


0.000048 


C: 


0.064 


0.0062 


0.0018 


0.0012 


N: 


0.070 


0.0049 


0.00053 


0.00043 


0: 


0.077 


0.0048 


0.00018 


0.00017 



9.5 Massive Fragments 

As it is mentioned above, spallation is characterized by observations of one 
heavy (in respect to the mass of initial target) nucleus, a small number of 
light fragments and several individual nucleons. Mass distributions of the 
heavy fragments, produced as result of proton induced spallation reaction on 
a chosen target vary with incident energy. In Figure 19.261 mass distributions 
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Table 9.5: Average values and corresponding standard deviations of mass 
of heavy nuclei remained after first and second stage of p+Pb reaction, at 
different incident energies (Tp) 



incident energy 


after first stage 


after second stage 


Tp = 0.5 GeV 


208.36±1.14 


188.21±7.17 


Tp = 1.0 GeV 


206.73±1.48 


186.65±7.31 


Tp = 2.0 GeV 


205.79±2.24 


188.33±9.70 


Tp = 3.0 GeV 


205.70±2.39 


189.53±9.98 


Tp = 4.0 GeV 


205.26±2.82 


188.68±10.24 


Tp = 5.0 GeV 


204.57±3.25 


187.70±10.58 


Tp = 7.0 GeV 


203.86±3.92 


186.89±11.28 


Tp = 9.0 GeV 


202.78±4.73 


185.04±11.83 



of heavy fragments formed in example p+Pb reactions, at several incident 
energies (calculated with the HSD + statistical evaporation model), set to- 
gether with mass distributions of corresponding residual nuclei remaining 
after first stage of the reactions (results of the HSD model only), are shown. 
It is seen, that both of the distributions (after first stage and after whole 
reaction) broaden with increase of incident energy. Additionally, the higher 
projectile energy, the lower average mass of nuclei after the first stage of the 
reaction. Distribution of mass of nuclei produced after the second stage of 
the reaction is much broader than mass distribution of the remnants of the 
first stage. The following changes of the distribution of mass of nuclei pro- 
duced after the second stage, as function of projectile energy are observed. 
At the lowest presented projectile energy, a maximum at low masses is seen 
and very small amount of nuclei with masses close to masses of nuclei after 
first stage of reaction (which are close to mass of the initial target). With 
increase of incident energy, the distribution broadens. The higher the inci- 
dent energy, the smaller is the maximum at lower masses. The lighter and 
lighter nuclei appears. Also more nuclei with masses close to mass of the 
initial nucleus are produced. Quantitative differences are presented in Table 

m 

In Figure 19.271 comparison of calculated mass and charge distributions 
of products of the proton induced reaction on Au target, at a few incident 
energies, with experimental data [70] and (HH!, is presented. Agreement be- 
tween the calculations and measurements is quite good. 
If one compares mass and charge distributions of products of proton induced 
reactions on a chosen target, at different projectile energies, significant dif- 
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ference is seen in mass range A > 25 and charge range Z > 10, as shown in 
Fig. 19.281 At the lower projectile energy, two almost separate distributions 
are visible: one (25 < A < 130, 10 < Z < 60, lower yield) corresponds to 
fission products, the second (130 < A < 200, 60 < Z < 80, higher yield) 
indicates spallation products (fission process is discussed in Sec. 18. 2p . At 
higher incident energy, the split of the distributions vanishes. 
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Figure 9.26: Mass distributions of heavy fragments remained after first and 
second stage of p+Pb reaction at few proton beam energies; results of the 
HSD (first stage) and PACE2 (second stage) model calculations (see also 
TabieESP 
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9.5 Massive Fragments 
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Figure 9.27: Mass (left) and charge (right) distributions of products of the 
second stage of p+Au reaction, at a few proton beam energies; results of the 
HSD+GEM model calculations (solid lines), compared to experimental data 
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Figure 9.28: Comparison of mass and charge distributions of products of the 
second stage of p+Au reaction emitted at different proton beam energies; 
results of the HSD+GEM model calculations (one observes that for higher 
incident energies of proton, the dips in the distributions become shallower) 
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Comparison of different models 
predictions on the first stage of 
proton induced reactions 

Apart from the HSD model, a few other microscopic models have been devel- 
oped in order to describe the fast stage of proton - nucleus reaction. All the 
models are based on different assumptions or simplifications (as described 
in Chapters [2] and [3]). It is interesting to compare their predictions. In 
this chapter, results of the HSD, INCL4 ([I7j) and QMD ([I9]) models are 
examined. 

The most important issue for comparison of results of various models of 
the first stage of proton - nucleus reaction are: particles production during 
the first stage of the reaction and properties of residual nuclei remaining 
after the stage. Predictions on pion production are particularly interesting 
for comparison. This is due to the fact, that calculations of the first stage of 
reaction are sufficient in order to receive realistic pion spectra, as all pions are 
produced only in the fast cascade of nucleon - nucleon collisions, where the 
available amount of four-momentum is large enough. One of the important 
features which can influence the magnitude of pion production and absorption 
is treatment of the nucleon-nucleus mean fleld (that corresponds to the main 
difference between the models). To check whether assumptions concerning 
the mean fleld in the model calculations affect signiflcantly multiplicity of 
produced pions, results of the HSD model are compared with the INCL4 and 
QMD models calculations and presented in Tables [TO. II and ri0.2[ for 2.5 GeV 
proton induced reactions on Ni and Au target, respectively. 

It is seen that the relative behavior of the pion multiplicities obtained 
from the different models is similar: for lighter target (Ni) the vr(+) : 7r(0) : 
7r(-) ratio is equal to 1 : 0.86 : 0.67 for QMD [82], 1 : 0.97 : 0.51 for INCL4 
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Table 10.1: Multiplicity of pions produced during p+^^Ni reaction at 2.5 GeV 
proton beam energy; results of the HSD, INCL4 fSTI and QMD [82^^ model 
calculations 





HSD 


INCL4 


QMD 


71+ 


0.51 


0.35 


0.42 




0.46 


0.34 


0.36 


7r~ 


0.38 


0.18 


0.28 



Table 10.2: Multiplicity of pions produced during p+^^^Au reaction at 2.5 
GeV proton beam energy; results of the HSD, INCL4 jSIJ and QMD jS^? 
model calculations 





HSD 


INCL4 


QMD 




0.49 


0.29 


0.34 




0.53 


0.35 


0.34 


7r~ 


0.46 


0.25 


0.32 



[8T] and 1 : 0.90 : 0.75 for the HSD, while for the heavy target (Au) 1 : 1 
: 0.94 (QMD [H^), 1 : 1.21 : 0.86 (INCL4 [H]) and 1 : 1.08 : 0.94 (HSD). 
However, absolute values of pion multiplicities obtained from the HSD model 
are always higher by factor ~1.5 than the values received from the INCL4 
and the QMD models, although results of the HSD model are closer to that 
of the QMD model. 

Let's compare now also an amount of kinetic energy carried out by the 
most abundantly produced particles during the first stage of the reaction, 
i.e. nucleons and pions. The average multiplicity and average kinetic energy 
of nucleons and pions emitted during the first stage of example 2.5 GeV 
proton induced reaction on Au target are displayed in Table [T073l The values 
collected in Table 110.31 indicate, that although the average multiplicity of 
emitted particles provided by the models are slightly different, the average 
amount of kinetic energy carried out by the particles is almost identical in 
all of the models. Consequently, properties of residual nuclei evaluated by 
the different models should be similar. In Figure 110. H comparison of the 
properties of residual nuclei, i.e. distributions of mass, charge, excitation 
energy and momentum in beam direction of the residual nucleus created in 
example p+Au reaction, at 2.5 GeV beam energy, are shown. The adequate 
average values are listed in Table 110. 4[ One observes that the HSD and 
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Table 10.3: Multiplicity, average values and standard deviations of total 
kinetic energy of nucleons and pions emitted during first stage of p+Au 
reaction, at 2.5 GeV proton beam energy; comparison of results of tie HSD, 
INCL4 mi and QMD |^ model calculations 





HSD 


INCL4 


QMD 




M 


<Ek> [GeV] 


M 


<Ek> [GeV] 


M 


< Ek > [GeV] 


n 


2.23 


0.78±1.27 


3.48 


0.80±1.46 


7.63 


0.76±1.83 


P 


2.20 


0.90±1.36 


2.68 


0.96±1.58 


5.63 


0.96±2.31 




0.49 


0.14±0.13 


0.25 


0.07±0.08 


0.32 


0.09±0.11 




0.53 


0.15±0.15 


0.35 


0.13±0.13 


0.34 


0.10±0.12 


7r+ 


0.46 


0.13±0.13 


0.29 


O.lliO.ll 


0.34 


0.11±0.12 



Table 10.4: Average values and standard deviations of properties of residual 
nuclei after first stage of p+Au reaction, at 2.5 GeV proton beam energy; 
comparison of results of the HSD, INCL4 fSlj and QMD [821 model calcula- 
tions 





HSD 


INCL4 


QMD 


< Ar > 


193.0±3.3 


191.8±4.03 


183.5±10.1 


<Zr> 


77.6±1.8 


77.3±1.9 


73.8±4.3 


<E* > [MeV] 


201.7±171.3 


253.6±226.3 


195.5±165.8 


< {Pz)r > [GeV/c] 


0.28±0.35 


0.44±0.51 


0.58±0.62 



INCL4 predictions are quite similar. 

Results of the models are sensitive to used parameters (in particular time 
duration of the first stage calculations, which is determined for each model 
individually, see [HI [19]). Therefore it is worthy to compare results obtained 
from each of the model and estimate their predictive power by confronta- 
tion with experimental data. In order to describe the nuclear reaction in 
unified way to verify predictions of the different models with experimental 
observations, the first stage models have been incorporated with one chosen 
statistical evaporation model applied for the second stage of reaction, i.e. 
GEM PU. Output of the HSD, INCL4 and QMD model, respectively (i.e. 
mass, charge, momentum and excitation energy of the residual nuclei) defines 
input for the GEM model. 

In order to verify which of the distributions of the properties displayed in the 



Chapter 1 0. Comparison of different models predictions 



149 




160 165 170 175 180 185 190 195 200 



80 8 




200 400 600 800 1000 -1 -0.5 0.5 1 1.5 2 2. 
E* (MeV) (Pz)r (GeV/c) 



Figure 10.1: Distributions of properties of tlie residual nuclei after first stage 
of p + Au reaction, at 2.5 GeV proton beam energy; comparison of results 
of the HSD (solid lines), INCL4 fSlf and QMD IS2S (scatter lines) model 
calculations 



Fig. 110.11 is the most reliable, calculated double differential kinetic energy 
spectra of ^He emitted at two different angles during p+Au reaction at 2.5 
GeV incident energy are compared with experimental data [7|, and shown 
in Fig. 110. 2[ Distributions of ^He are the best choice for such a comparison 
due to the fact, that in the ^He channel evaporation is dominant. As it is 
shown in Ref.[7], about 90% of the "^He spectrum consists of the particles 
produced during the evaporation stage. One can see, that the models give 
similar results, although the HSD and the QMD model results are somehow 
closer to experimental data. 
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In Figure 110.3^ comparison of calculated example fragment mass and charge 
distributions compared with experimental data [TU] and [20] sue displayed. 
It is seen, that also in this case the models give similar results. In particular, 
results of the HSD and the INCL4 models are almost identical. If looking at 
the QMD model calculations, plotted for reaction p+Au at 1.0 GeV beam 
energy, in the emitted fragment mass number and charge range correspond- 
ing to the fission peak (i.e. 50 < A < 100 and 20 < Z < 50), they are lower 
than calculations of the other models, but closest to the experimental data. 
However, outside the range (i.e. A < 50, A > 100 and Z < 20, Z > 50) 
results of all considered here models are very similar. 
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Figure 10.2: Double differential distributions of'^He emitted during p+Au re- 
action, at 2.5 GeV incident energy; results of the HSD+GEM, INCL4+GEM 
(I81]j) and QMD+GEM (IS'Jj) model calculations compared to experimental 
data measured by PISA collaboration |^ 



Based on the comparisons, one can conclude that results of the HSD 
model do not differ significantly from results of the other models. The addi- 
tional verification by confrontation of the calculations with experimental data 
indicates, that results of the HSD and the QMD model are particularly close 
to the experimental data. As results of the HSD model are not distinctly 
different from that of the QMD model (which is actually the most advanced 
one), a big advantage of the HSD model, connected with respectively short 
time of calculations should be mentioned. The average time needed to cal- 
culate a reaction with the QMD model is of the order of a few days, while 
time needed to calculate the adequate reaction with the HSD model - only 
few hours. 
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Figure 10.3: Comparison of mass (left) and charge (right) distributions of 
products of the second stage of p-hAu reaction emitted at different pro- 
ton beam energies; results of the HSD+GEM, INCL4+GEM (l81l) and 
QMD+GEM (for Tp = 1.0 GeV only, 182]) model calculations compared 
to experimental data and fSUI 
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Summary and Conclusions 



Global information of properties of proton - nucleus reactions (i.e. particles 
production, properties of residual nuclei after first stage of the reactions) has 
been presented in wide range of projectile energy and mass of target nuclei, 
with the Hadron String Dynamics (HSD) [381 EH] model for the first stage of 
the reaction, incorporated with statistical evaporation model (GEM [2T] and 
PACE2 [SU]) for the second stage. It has been confronted with experimental 
data and results of other models. 

The version of the HSD code ^38j has been developed in order to describe 
proton induced spallation reactions in wide energy range in the following 
aspects. First of all, to calculate properties of residual nuclei after first stage 
of the reaction, defining an input for models of the second stage calculations. 
The code has been developed for description of pion production in proton 
- nucleus reactions. Additionally, a version of HSD code that allows for 
calculations of pion induced reactions has been prepared. 
Stopping time for the first stage model calculations has been established 
based on behaviors of excitation energy, longitudinal momentum and angular 
momentum as a function of time. These are decisive for determination of the 
duration time of the first stage of proton - nucleus collisions, due to the 
fact that in the dependences two phases can be distinguished: a phase of 
rapid variations, followed by a phase of much slower variations, indicating 
thermodynamical equilibrium of the nucleus. Properties of residual nuclei 
are evaluated by exploring conservation of energy, mass, charge, momentum 
and angular momentum. 

It has been shown, that distributions of properties of residual nuclei (i.e. 
mass, charge, excitation energy, momentum, angular momentum) differ sig- 
nificantly with mass of target, but behave similarly for varied projectile en- 
ergies. Multiplicity of emitted pions increases with incident energy and with 
target mass number. Contrary, e.g. average kinetic energy of other individ- 
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ual emitted particles (i.e. nucleons, pions) increases with kinetic energy of 
incoming proton, but decreases with target mass number. Higher nucleon 
and pion momenta and their increasingly forward peaked distributions are 
associated with increasing incident energy. That is connected with increas- 
ingly forward focused nature of nucleon - nucleon collisions with increase of 
incident energy; the total nucleon - nucleon scattering cross section remains 
approximately constant, but ratio of inelastic to elastic collisions increases 
rapidly with impact energy |45j . 

It has been shown, that the biggest amount of initial energy deposited 
in nucleus by incoming proton is carried out during the first stage of the 
reaction, i.e. about 85 %, and only about 15 % of the initial energy is 
carried away during the second stage of the reaction. Most of the initial 
energy is carried out by nucleons and pions. Neutrons are emitted mainly 
in the second stage of reaction, whereas protons - mainly during the fast 
cascade stage. Pions are emitted only in the first stage of reaction, where 
the locally available amount of four-momentum is large enough. In case of 
light composites, only low energy part of distribution, corresponding to a 
result of evaporation stage, is described by the HSD plus evaporation model. 
E.g. for deuterons, the absence of high energy part of modelled distributions 
indicates a lack of deuterons produced during first stage of the reaction. 
Description of the whole distributions necessitates future implementation 
of some additional mechanisms into the first stage model, e.g. coalescence 
processes. It is observed already in neutron spectra that these additional 
mechanisms have to depend strongly on size of the system: for heavier targets 
the neutron spectra in some region are visibly underestimated by model 
calculations, see Fig. 19.101 

Very interesting conclusion of the work is a prevision for a negligible 
probability for fragmentation in proton induced reactions, at considered here 
range of incident energy. This statement is based on results of the HSD 
model calculations of excitation energy per nucleon of residual nuclei and 
time evolution of nucleon density. Results of the HSD model show, that a 
value of excitation energy greater than 5 MeV/N occurs very rarely, in only 
few percent of cases (studies of fragmentation have found that it occurs at a 
value of excitation energy per nucleon of a residual nuclei at least equal to 
about 5 MeV/N [Ml [55] ). Fragmentation could be more probable if a part of 
nuclei takes part in excitation, i.e. the excitation energy should be calculated 
only per some amount of the active nucleons; only this part of nucleus would 
undergo fragmentation. However modelling of nucleon density evolution with 
the reaction time shows, that proton induced reactions are very low - invasive 
processes. Incoming proton causes insignificant modifications of the nuclear 
density. Energy is rather uniformly distributed among the nucleons of resid- 
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ual nuclei. A situation, when two differently excited parts of nucleus appear, 
occurs only in about 1% of proton - nucleus reaction cases and rather in 
peripheral collisions, whereas experimentalists associate fragmentation with 
central collisions. 

It has been shown that in proton induced reaction, in the considered here 
incident energy range, spallation is evidently dominant process. Moreover, 
the second stage of spallation reaction is a competition between evaporation 
and fission. 

Results of the HSD model calculations have been compared with available 
experimental data and with adequate results of other microscopic models 
(i.e. INCL4 [T7j and QMD [19]). The comparisons model - model show 
that results of the HSD model do not differ significantly from results of the 
other models calculations. Results of the HSD model are also consistent 
with various experimental data. Qualitatively, the data are well described 
by the calculations. Quantitatively, the maximal noticed discrepancy is of 
factor 2. It concerns pion spectra, which are the worst described cases, i.e. 
overestimated by a factor of 2. That is because, despite description of NN 
interactions is well established, details of pions propagation in nuclear matter 
are very delicate. Legth of free path for pions in nuclear matter of standard 
density is 1 - 1.5 fm. It means, pions that arrive to detectors stem from 
surface of nucleus. In models like HSD it is just the surface that is described 
only approximately. Pions production due to proton - nucleus collisions is 
a complicated process, since they are not emitted directly, but result from 
a succession of productions and absorptions inside a nucleus, mediated by 
Delta resonance. 

The comparisons indicates that assumptions and simplifications employed 
in the HSD model are correct. 

In order to obtain the possible best description of proton induced reac- 
tions, trials of various feasible modifications of the HSD model parameters 
has been investigated. However, the possible modifications, as in case of each 
of semi-classical models are limited. A necessity of the exact fully quantum 
model appears. Unfortunately, construction of such model of many body 
strongly interacting system is at present practically impossible. 

Summarizing, the global properties of proton induced reactions have been 
presented in wide range of mass of target and incident energy. It has been 
shown, that these are quite non-invasive reactions, where spallation is the 
dominant process. Properties of residual nuclei after the first stage of the 
reaction are weakly dependent on incident energy, but strongly dependent 
on mass of target. High momenta of emitted particles are associated with 
high incident energies. The second stage of the reaction is a competition of 
evaporation and fission. The agreement of results of the HSD calculations 
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with experimental data and with results of the other microscopic models 
calculations indicates that proper assumptions have been employed in the 
HSD model. As discussed in Sec. 15.11 results obtained from the HSD model 
vary smoothly as function of proton incident energy and mass of target, so 
interpolation of results is quite feasible. 
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